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ABSTRACT 


A readiness-based sparing (RBS) model for the repair and replenishment of repairable items is needed 
by the Navy which considers the aggregate inventory of repaired units and new ones. This thesis presents 
progress in the development of such a model. In contrast to other such current repairables models in the 
literature, it also allows for both batch repair and procurement. A theoretical model had been developed 
earlier at the Naval Postgraduate School for the probability distribution of inventory position for such a 
model. However, no theoretical model has yet been developed for the probability distribution of net 
inventory because the real-world inventory management of repairables is quite complex. Therefore, a 
simulation model was developed of the Navy's repairables management process to explore the nature of that 
distribution as a function of relevant system parameters. It was then run for a range of values of a subset 
of those parameters. The net inventory distribution appears to be normally distributed with its mean and 
variance being a linear function of the product of carcass return rate and repair survival rate. The theoretical 
distribution for inventory position was not only validated, it was found to be quite robust. Further analyses, 


however, are required before the effects of all relevant parameters are well understood. 
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EXECUTIVE SUMMARY 
A. PURPOSE 

The purpose of this thesis is to develop a simulation 
model which accurately describes the Navy’s secondary item 
repairables process, and then use the simulated data to 
study the aggregate net inventory and inventory position 
distributions, compare the resulting inventory position 
distribution to a proposed theoretical one, analyze the effect 
of three parameters on the mean and variance of the net 
inventory, and define safety level. 

The main reason for developing a mathematical model of an 
inventory system is to use it as an aid in developing an 
operating doctrine. Because the mathematics in the 
repairables model are so complicated, to accurately model it 
analytically is virtually impossible except in the limiting 
case where the repair induction quantity is one. However, 
simulation can be used to develop an approximate analytical 
model. 

B. BACKGROUND 

Both Navy ICP’s use the Uniform Inventory Control 
Program (UICP) system for determining the replenishment of 
consumable and repairable secondary items. The formulas used 
in the UICP for determining the order quantity and reorder 


point for the repairables replenishment system were developed 











in the late 1960's. However, no rigorous mathematical 


development supported them. Rather, they were direct 
extrapolations of the consumable model formulas. The result 
was two sets of formulas, one for the procurement process and 
one for the repair process. In particular, there were two 
separate sets of risk equations, safety levels, order 
quantities, and reorder levels. The only "improvement" to the 
repair model occurred in 1984 when the Integrated Repairables 
Model was implemented. This model contained only one risk 
formula and applied the same safety level computation to both 
processes, but continued to model the repairable inventory 
management process as a two-part process. 

C. RESULTS 

Two major results of this study have emerged. The first 
is that the theoretical distribution for inventory position 
was not only validated, but found to be quite robust. The 
second is that the net inventory distribution appears to be 
normally distributed with its mean and variance being a linear 
function of the product of the carcass return rate and the 
repair survival rate. 

D. RECOMMENDATIONS. 

More sensitivity analyses of the inventory position and 
net inventory probability distributions must be accomplished. 
Clearly, the effect of the low values of the carcass return 
rate on the net inventory distribution parameters needs to be 


better understood. In addition, major sensitivity analyses of 


xii 








these distributions must include the effects caused by varying 
demand rates and distributions, and by varying turn-around 


time rates and distributions. It is also recommended that the 


Naval Postgraduate School and the Naval Supply Systems Command 


share in the research effort so that it is accomplished in an 


expeditious manner. 
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I. INTRODUCTION 


A. BACKGROUND 

As of 30 September 1992 the Navy’s two Inventory Control 
Points (ICP’s), the Aviation Supply Office in Philadelphia, PA 
and the Ships Parts Control Center in Mechanicsburg, PA, 
managed over 580,000 secondary support items which were broken 
down into two major classifications, consumables and 
repairables [Ref. 1, 2). Consumable items are ones which are 
immediately disposed of at the time of failure. Repairables 
are items which are sent to a repair activity upon failure and 
are typically returned to "new" condition at a cost which is 
no more than that of a new item and in a time interval 
substantially shorter than a procurement lead time. 

As a result of the Defense Management Review Decision 901, 
"Reducing Supply System Cost/Navy Inventory Reduction Plan", 
the Navy is in the process of transferring almost all of its 
consumable items to the Defense Logistics Agency leaving the 
ICP’s with some 207,000 repairables to manage. To maintain 
its existence, the ICP’s need to better manage this material, 
i.e, each must become a repairables "Center of Excellence". 

Both Navy ICP’s use the Uniform Inventory Control Program 


(UICP) system for determining the replenishment of consumable 











and repairable secondary items. The formulas used in the UICP 
for determining the order quantity and reorder point for 
consumables were developed in the mid-1960’s based on the 
research of Professors G. Hadley and T. M. Whitin. [Ref 3] 
The formulas for the repairables replenishment system were 
developed in the late 1960's. However, no rigorous 
mathematical development supported them. Rather, they were 
direct extrapolations of the consumable model formulas. The 
formulas for procurement and repair order quantities and the 
associated reorder points were based on the modelers’ 
intuition from modeling consumables. The result was two sets 
of formulas, one for the procurement process and one for the 
repair process. In particular, there were two separate sets 
of risk equations, safety levels, order quantities, and 
reorder levels. The only "improvement" to the repair model 
occurred in 1984 when the Integrated Repairables Model was 
implemented. This model contained only one risk formula and 
applied the same safety level computation to both processes, 
but continued to model the repairable inventory management 
process as a two-part process. 

The development of an analytically-based repairable item 
inventory system is an extremely complex process. As will be 
seen in the discussion of the development of the Navy’s 
current model, when the formulas did not provide values which 


the inventory managers could believe in, more assumptions were 


made, and more approximations were developed. However, the 











problems were really rooted in the basis of the model being 


incorrect, and subsequent “patches" did not help much. 


B. OBJECTIVE 

The main reason for developing a mathematical model of an 
inventory system is to use it as an aid in developing an 
operating doctrine. Because the mathematics in the 
repairables model are so complicated, to accurately model it 
analytically is virtually impossible except in the limiting 
case where the repair induction quantity is one. However, 
simulation can be used to develop an approximate analytical 
model. Indeed, Hadley and Whitin suggested that in very 
complex cases one may have to use a computer and simulation to 
study a small number of different operating doctrines and 
choose one which best suits the objective. (Ref. 3] 
Therefore, the objectives of this thesis are to: 


1. Develop a simulation model which accurately describes the 
Navy’s secondary item repairables process. 


2. Use the simulated data to study the aggregate net 
inventory and inventory position distributions. 


3. Compare the resulting inventory position distribution to 
a proposed theoretical one. 


4. Analyze the effect of the Carcass Return Rate (CRR), the 
Repair Survival Rate (RSR), and the Repair Processing 
time (REP) on the mean and variance of the net inventory. 


5. Define safety level. 


C. SCOPE 

The intent of this thesis is to take the first step in the 
analysis of the Navy’s aggregate repairables inventory system. 
The simulation analysis focusses on the _ probability 
distributions for inventory position and net inventory. 
However, the scope of the simulation studies includes only 
three stochastic parameters. One stochastic parameter is the 
quarterly demand, and it is assumed to be generated by a 
Poisson process. The two other stochastic parameters are the 
Carcass Return Rate (CRR) and the Repair Survival Rate (RSR) 
which are probabilities and are generated by a Bernoulli 
process. All other parameters are assumed to be deterministic 
because the models described in Reference (1) are founded on 
the same assumptions and are serving as a basis for 


comparison. 


D. PREVIEW 

Chapter II is an indepth discussion of the development of 
the inventory system formulas for the Navy’s Consumable Model 
and the Navy’s Integrated Repairables Model. Chapter III 
introduces the inventory position probability mass function 
formulas for the proposed theoretical model, and then 
describes in detail the simulation model, the way data was 
generated, and how the data was manipulated. Chapter IV 


contains the analysis of the results of the simulations. 





Chapter V summarizes the thesis research and findings, offers 


conclusions, and makes recommendations for further analyses. 








II. DEVELOPMENT OF UICP INVENTORY MODELS 

Inventory problems are as old as history itself, but it 
has only been since the beginning of this century that 
mathematics have been used to manage inventory systems. These 
attempts focused on what will be called "the operating 
doctrine"; namely, how much to buy (order quantity) and when. 

The earliest of these mathematical “models" was developed 
by Ford Harris of Westinghouse Corporation and the order 
quantity he derived was called the "simple lot size formula" 
(Ref. 3]. This formula was subsequently used after World War 
I by R. H. Wilson and ‘is often referred to as the "Wilson 
Formula" because he used it as an integral part of the 
inventory control scheme he sold to many organizations. The 
model assumed demand to be known and constant, and no stock- 
outs were allowed. 

Following World War II, inventory models which considered 
demand to be a random variable were developed [Ref. 4]. Among 
those who developed these models were G. Hadley of the 
University of Chicago and T. M. Whitin of the University of 
California, Berkeley. They received support from the Navy’s 
Bureau Of Supplies and Accounting, the predecessor to the Navy 
Supply Systems Command (NAVSUP), for much of their work. In 
1963 this research was consolidated in their definitive text, 


"Analysis of Inventory Systems", which was published by 


Prentice Hall, Inc., Englewood Cliffs, New Jersey. Their 
research formed the cornerstone of the Navy's wholesale 
inventory models, [Ref. 5], and in particular, the Uniform 
Inventory Control Program (UICP) consumable item inventory 
model, (Ref. 6]. The next section presents the derivation of 
this model and, because it forms the basis of the Navy’s 
current wholesale inventory system, it will serve as a 
reference for the models described later for managing 


repairable items. 


A. THE NAVY’S UICP CONSUMABLE MODEL 
Basic to the Navy’s UICP consumable model are the 
definitions of Inventory Position (IP), Net Inventory (NI), 


and the average annual Total Variable Costs (TVC) where, 


IP = On Hand + On Order - Backorders 
=O/H + 0/0 - B/O; a) 
NI = On Hand - Backorders 


O/H - B/O; (2) 


and 


TVC = Average Annual Order Costs + 
Average Annual Holding Costs + (3) 
Average Annual Backorder Costs. 











are: 


10. 


il. 


The major assumptions of the Navy’s consumable model 


. There exists a steady state environment. This means that 


demand, while being a random variable, has a fixed mean, 
variance, and probability distribution which do not 
change over time. 


. There is a transaction item reporting system; whenever 


a demand occurs, it is immediately recorded. The result 
is a continuous review of demand. 


. Units are demanded one at a time. This assumption is 


used to develop the order costs and holding costs terms 
of the model (they are taken directly from Reference 3). 
The UICP model violates this assumption in the backorder 
costs term and allows a requisition size of any number of 
units. 


. Ordering occurs when the inventory position reaches a 


value called the reorder level or reorder point. That 
value is non-negative. 


. The amount to be procured is called the order quantity 


and has a constant value. Procurement of the order 
quantity is not constrained by budget restrictions. 


. The cost of placing an order is constant and independent 


of the order quantity. Procurement lead time can be 
either deterministic or a random variable. In this 
section it will be assumed to be known and constant. 


Inventory holding costs are proportional to the unit 
cost. 


. Demands are either filled or backordered. 


. The cost of a backorder is quantifiable. 


A factor, known as Military Essentiality, is 
quantifiable. This factor attempts to measure an item’s 
worth in terms of operational readiness. An item such as 
an air/search radar would have a large factor because 
loss of one would severely cripple a ship’s mission, 
whereas, night-vision goggles would have a smaller factor 
because loss of one would not cripple a ship’s mission. 


An entire order quantity is received at the same time; 
i.e., it is not split up into two or more subsets 
arriving at different times. 








The objective of the model is to determine non-negative 


values of the order quantity, Q, and the reorder point, R, 


which minimize equation (3). 


B. DERIVATION OF THE TVC FORMULA 
1. Average Annual Order Costs 
The average annual ordering costs are simply the 
product of the cost of placing an order and the expected 
number of orders placed per year. The cost of placing an 
order is denoted as A. The expected number of orders placed 
per year is the forecasted average annual demand, 4D, (where 
D is the forecasted quarterly demand), divided by Q. Thus, the 


average annual ordering costs can be written as equation (4), 


Average Annual Order Costs = A : (4) 


2. Average Annual Holding Costs 
The average annual holding costs are computed as the 
product of the average on-hand inventory and the annual 
holding cost rate per unit, IC, where C is the unit purchase 
cost and I is a percentage established by the Fleet Macerial 
Support Office and is composed of four estimated factors, a 
"value of money" rate of 10%, an obsolescence rate of 10%, a 


shrinkage/pilferage rate of 2% anda "handling" rate of 1%, 








for a total of 23%. The average on-hand quantity, E(O/H), is 


found by rewriting equation (1) as 


E(O/H) = E(IP) - E(0/0) + E(B/O). (5) 


The model assumes that every time a demand occurs, the 
IP is reduced by one unit. At the instant when the inventory 
position reaches the reorder point, R, an order of size Q is 
generated and the inventory position, IP, increases to R+Q 
because Q has been ordered. Hadley and Whiten [Ref. 1:p 183) 
have shown that, under the assumption of Poisson demand, the 


probability of IP being R + x is 


p(IP=R+x) = 4 forlsxsQ. (6) 


Thus, the expected IP at any point in time is 


E(rp) = Q*{R*R+2 


(7) 





The expected on-order quantity can be found in the 
following way. The time between generations of an order is 
known as the order cycle. Let M be the expected number of 
order cycles in a procurement lead time, L. Then, 


L = (M (order cycle). (8) 


Multiplying both sides by the forecasted quarterly demand, D, 


results in 


10 

















(D) (Z) = (M (order cycle) (D). (9) 


If the order cycle begins with the generation of an order, 
then another order will be generated when the demand reaches 


Q. Therefore, 


DL = MO, (10) 
where MQ is the expected on-order quantity. The result is 
that the expected on-order quantity is equivalent to the 
expected lead time demand. 

Substituting equations (7) and (10) into (5) and 
denoting E(B/O) as B(Q,R), the resulting equation for E(0/H) 
is 


E(O/H) = ent + R- DL + B(O,R). (11) 


The average annual holding costs are, therefore, 
Average Annual Holding Costs = IC (2+ +R- DL + B(0.R)). 
(12) 


3. Average Annual Backorder Costs 
The UICP determines the cost of backorders based on 
the expected time-weighted number of requisitions short per 
year. This is done by determining the expected number of 
backorders on the books at any time t, and then dividing that 


value by the average requisition size, S. 


11 


Ee 


To find the expected number of backorders at any time 
t, let R + x be the inventory position at the beginning of a 
procurement lead time preceding t (i.e., at time t - L) and 
let y equal the number of backorders occurring during that 
lead time. To incur y backorders, the demand over the lead 
time must be R + x + y since any quantity on order prior to 


the beginning of a lead time will have been received by the 





end of the lead time, and any order placed after the beginning 
of an order cycle will be received after the end of the lead 
time (as noted above, lead time is assumed to be known and 
constant). Let the Poisson probability mass function of lead 
time demand u be denoted by p(u;DL) with the corresponding 
cumulative distribution function, P(u;DL). The combined 
probability that R + x + y demands occur during a lead time 
and the inventory position was R + x at the sta:t of the lead 


time is the product 


P(R+x+y;DL) p(IP=R+x), (13) 
and, because of equation (6), (13) reduces to 
3 p(R +x +y;DL). (14) 


Equation (14) is valid for all x such that 1 s x s Q, and all 
y such that 0 s y s o. 
To find the probability f(y) of y backorders for all 


possible IP values, equation (14) must be summed over all the 


possible values of x: 





Q 


yy 2ptR + x + y;DL) 


f(y) 
yey 2 


3 (PUR + 0+ y:DL) - P(R + y;DL)}; (15) 


and the expected number of backorders at any time t, B(Q,R), 


is the expected value of y: 


B(O,R) = > yf(y) 


ys0 


YD y(4) (P(R + 0+ y;DL) - P(R + y;DL)] 
y=0 o 


10 yip(R+Q+y:DL) - P(R+y;DL)) . — ¢46) 
y=0 


If we let u=R + y, then (16) can be rewritten as 


B(O,R) = 2D (u- R)[P(u+ Q;DL) - Plu;DL)] . (47) 


OF 


As noted in Reference (3), B(Q,R) is also the average 
unit-years of backordered units per year. To get the average 
time-weighted annual costs of requisitions short, the UICP 
model first divides equation (17) by the average requisition 
size, S, and then multiplies the result by the product of a 
shortage cost per requisition, >A», and the military 


essentiality factor, E. The result is 


Average Annual B/O Costs = (AZ}5(0,x) : (18) 
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4. Average Annual Total Variable Costs 
Substituting equations (4), (12) and (18) into 
equation (3) yields the following formula for the average 


annual Total Variable Costs: 


4aDA 
Q 


+ 1c( 24 + R- DL + B(O,R)) + 4EB(0,R) (19) 


5. Determining the Optimal Values of Q and R 
To simplify discussion of the optimization process and 
the development of the associated formula for Q and R, lead 
time demand will be assumed to be large enough that a Normal 
distribution can be used to approximate the Poisson. This 
will allow us the use of calculus to determine Q and R. In 
the process, F(u;DL) will denote the cumulative distribution 


function for the Normal distribution. In this continuous 





case, 


B(Q,R) = Zf(u-R) [F(u + Q;DL) - F(u;DL)]du. (20) 
R 


Taking the partial derivative of (19) with respect to 


Q results in 


arve , -4DA , IC AE) Oo 
SRO ot tet EC + AF) Zp BlORd. (21) 





Unfortunately, 210.8) is impossible to write out as a 


simple general function of Q because Q is contained in 
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F(u + Q;DL). As a consequence, this term is ignored in UICP 
and the remainder of the partial derivative is set to zero to 


get 


+— =09Q, 


o?z = 2 





which, when solved for Q°, gives: 


or =| 228. (23) 


After Q° is computed, the UICP constrains the order 
quantity to at least one quarter’s worth of demand and no more 
than six quarters’ demand. It is further constrained to avoid 
any possible deterioration due to shelf-life. 

The next step is to determine the optimal reorder 
point. This is done by first determining the optimal 
probability of stockout (the UICP calls it "RISK") and 
comparing it to the acceptable range constrained in the UICP. 
Once the RISK value has been determined, the reorder point can 
be computed based on the probability distribution of the lead 
time demand. 

The RISK formula is developed by taking the partial 
derivative of equation (19) with respect to R and setting the 


result equal to zero, 


arve _ a AE| O 7 
oe tice 1c (3,810, R)) + AS BRE (OR) 0. (24) 
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Equation (24) reduces to 


IC«- (zc * AE BIO.) . (25) 


The partial derivative 2B(0.R) is obtained using 


Leibnitz’s rule for differentiation of an integral. The 
result is 
© B(0,R) = -Af tru + O;DL) - F(u;DL)) du. (26) 
OR oh ’ ’ 


Substituting (26) into (25) results in 
rc = (rc + F\afre + Q;DL) - F(u;DL)] du. (27) 


Multiplying both sides by the product, (S) (Q), and dividing by 


(SIC + AE) results in 


fra + 0;DL) - F(u;DL)]du = 
R 


ICOS 
SIC + AE (28) 
Equation (28) is impossible to solve analytically for 
R since the left-hand side cannot be reduced to a simple 
function of R. An alternative approach is to use an 
approximation technique, described in Reference (7), to obtain 
the RISK formula. Using some limiting arguments, it can be 


shown that 
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- F(R; — SIC __ = ; 
O(1 - F(R;DL)) > = 2 Ol - F(R + Q;DL)), (29) 


which reduces to 


rs , SIC es ; 
(1 - F(R;DL)) 2 Sic eur? (1 - F(R + 0;DL)) (30) 


Then, instead of solving for the smallest R that satisfies 


equation (30), UICP merely uses the equation 





to solve for R. This expression describes the probability of 
being out of stock during a lead time, or RISK. The final step 


is to let 


(32) 


¢,) 
1] 
=lb 


where W is the average quarterly frequency of requisitions. 


Substituting (32) into (31) gives: 


DIC (33) 


RISK =1 - F(R; DL) = Dic + wEW: 


UICP then places both a lower and upper bound on the value 
that RISK can take. The bounds depend upon how active the 
item is; for instance, at the Navy’s Ships Parts Control 
Center, most of the DLR items they manage have bounds of 0.15 
and 0.4, respectively. 

The <eorder level, R, is obtained by solving equation 


(33), constraining RISK with an upper and lower bound, and 
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then determining the probability distribution of the lead time 
demand. The distribution assumed is determined by two 
methods. One method involves the category of the material. 
UICP divides all wholesale material into categories called 
Marks. A Mark 0 item is one in which the yearly demand is 
very small, usually less than one. The lead time demand for 
all Mark 0 items is assumed to be Poisson distributed. For 
all items in the remaining categories, the ICP’s have 
established a value known as the probability break-point. If 
the lead time demand is greater than or equal to this value, 
it is assumed to be Normal. Otherwise, it is assumed to be 
Negative Binomial. 

There are three methods in calculating R. When the 
lead time demand is assumed to be Normal, the reorder point 


becomes the lead time demand plus the safety level: 


R= DL + Zo, (38) 


where o is the standard deviation of demand during lead time, 
and z is the standard Normal deviate associated with the value 
for RISK. When the lead time demand is assumed to be 
Poisson, R is the smallest integer which satisfies the 


following inequality: 
F(R) 21 - RISK. (35) 


This same method is employed when lead time demand 


distribution is assumed to be Negative Binomial. 
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C. THE NAVY’S DEPOT LEVEL REPAIRABLES MODEL 

As noted, the Navy’s consumable inventory model conforms 
to the requirements of Reference (5). However, there is no 
such guidance for the procurement or repair of Depot Level 
Repairables (DLR) because the group that developed Reference 
(5) realized there would be little hope for consensus on a 
repairables document. That is because each service had 
developed its own model with no two models being exactly 
alike. 

This section will focus on the Navy’s model. The 
objective of the Navy’s model is to determine how much to buy, 
when to buy, how much to repair, and when to repair so that 
the average annual total variable costs are minimized. 

The Navy’s DLR model is based on the consumable model with 
a couple of major differences. First, the UICP views the DLR 
system as two separate systems in the modeling process; one 
for the procurement of new material and the other for the 
repair of the not-ready-for-issue (NRFI) but repairable 
carcasses. Thus, the system receives ready-for-issue (RFI) 
material from two sources, repair and procurement, where the 
repair rate G (called the "regeneration" rate in UICP) is 
forecasted and the attrition rate is (D - G). The 


regeneration rate is computed as 
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G = (CRR) (RSR) (D), (36) 


where CRR is the estimated carcass return rate, RSR is the 
forecasted carcass repair survival rate, and D is the 
forecasted quarterly demand. 


Second, after a certain number of attritions occur, a 


procurement is needed to buy replacements for the lost units. 


The lead time attrition demand for the “procurement system" is 





assumed to be composed of a weighted average lead time L, and 
the quarterly demand rate D. It is referred to as the 
procurement problem variable (Z). From Reference (6), L, is 


given by equation (37): 


-G G 
t+ (25°) + Bh ae 
where T, is the repair turn-around time. The ratio ip 6) 7 S) 


is the fraction of the time that a demand for an RFI unit has 
no associated repairable carcass being returned to a depot. 
Conversely, the ratio < represents the fraction of the 
time that a demand is accompanied by a repairable carcass. 
If a repairable carcass is not returned, a new unit will 
have to be purchased. That unit will take a procurement lead 
time to be received by the system. Thus, L is associated with 
that fraction of demands having no repairable carcasses. When 


a carcass is returned and is repairable, it will take a time 
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T,, called the repair turn-around time, for it to be processed 


through repair. 


The equation for the Procurement Problem Variable, Z, is: 


Z = DL, = (D- G)(L) + Gh. 138) 


The lead time demand for the “repair system", Z, is 


denoted by the product of D and the repair turn-around time, 
Z, = DT;. (39) 


Repair turn-around time is measured from the time a failed 
carcass is inducted into the repair system (it changes from 
condition code F to M) until it is successfully repaired and 
returned to RFI condition (from M to A condition). 

Since the UICP views the DLR system as two separate 
systems in the modeling process (one for the procurement of 
new material and the other for the repair of the NRFI but 
repairable carcasses), two inventory positions are used, [Ref. 
6]. The inventory position for the procurement problem, IP,, 


(ignoring any planned requirements and war reserves) is: 


IP, = (On hand RFI material) 

(RFI on order through procurement) 

(RSR) (On hand NRFI carcasses) (40) 
(Backorders) 

(RSR) (Unfilled NRFI requirements for repair). 


r+ 
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The inventory position of the repair system, IP,, is: 


IP, = (On hand RFI material) 


+ (On order procured material 

Scheduled to be received within 

a depot level repair turn-around time) (41) 
+ (RSR) (Inducted material) 
- (Backorders) . 


1. Procurement Problem 
The average annual total variable cost equation 
associated with procurement was adapted from the consumable 


model by incorporating the attrition demand (D-G) and the 


Procurement Problem Variable, Z, into equation (19). The 
result is 
4(D-G)A AE(4(D-G 
TVC = Ate la + IC (2 +R-Z+ B,] + B(te Oe. (42) 


where B, is the expected number of backorders at the end of a 


lead time. The formula for B, is given by Reference (6) as 


B, = [(u- R)£(u;DL) du (43) 
R 


because, when the model was developed, there was no 
requirement for the expected backorders to be time-weighted. 


Taking the partial derivative of (42) with respect to Q gives 


OTvc . -4(D-G)A , IC , ~4AE(D - G)B 
arom - eer aes 3. ae (44) 


When equation (44) is set equal to zero, the result is: 
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oe | 8(D- G)A , BAE(D - G)B, (45) 
57a, ICS 


Unfortunately, B, is a function of the procurement problem 


reorder point R. Thus, to solve (45) would require an 
iterative procedure which includes determining optimal R as a 
function of Q. Because of this, the UICP approach is to 


ignore the second term and use 


ors | 8(D- G)A (46) 
Ic ° 


The order quantity constraints mentioned before for the 
consumable model are .then applied after replacing D by 
(D - G). 
2. Repair Induction Quantity 
The average annual total variable costs associated 
with the repair system were only "surmised" as indicated in 
Reference (6). The formula used is: 


rve, = AIMIn(D. GI g, « ro[% +2, -2,+2,] 


4A E[Min(D, G) 
* SREVEDIE SI a, (47) 
2 
where C, is the unit cost to repair an item, A, is the cost to 
prepare a repair order, Q, is the repair quantity, R, is the 
repair reorder point, B, is the expected number of units 


backordered during T,, B, is the expected number of 
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requisitions backordered during T,, and Min(D,G) is the 


smaller of the forecasted quarterly demand or the quarterly 
regeneration rate. The equations for B, and B, are not 
presented here because neither is used in the integrated 
repairable model. 

The reason presented for using Min(D,G) is that it is 
possible that forecasted regenerations exceed forecasted 
demand. (Examples include a phase-out of a weapon system 
where there exists turn-ins but little or no demand, or where 
an extreme decrease in one quarter’s demand coincides with the 
receipt of a large number of carcasses from demands of the 
previous quarter.) 

Following the same approach used in the consumable and DLR 


procurement models, the optimal repair order quantity is 


os 8 [Min (D,G) ]A, (48) 
$ zc; , 


and is constrained with the same limits. 


approximated as 


3. Integrated RISK Formula 
In the past, the procedure for computing reorder 
levels separately for the DLR procurement and the DLR repair 
quantities resulted in the procurement levels for many items 
not being sufficient to provide carcasses so that the repair 


quantity Q, could be inducted. [Ref. 6] This was due, in 
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part, to the computation of two different safety levels and, 
in certain cases, the repair safety level exceeding the 
procurement safety level. This motivated the integration of 
the computation of the safety levels for the procurement and 
repair levels. The first step of the process was to compute 
a weighted average unit/repair cost, C,, from C and C, in the 
same way as L, was developed, and then substitute C, for C in 


equation (33). The result is 


DIC, 


RISK=1 - F(R; DL) = DIC, + AEW’ (49) 


As with the consumable case, the RISK is constrained by UICP 
set minimum and maximum values. 

The procurement reorder point, R, is computed in the 
same way it was for the consumable model. The associated 


safety level follows from the fact that: 


R=2Z+ safety level. (50) 

Once the safety level is determined, the repair reorder point, 
R,, is calculated using the same safety level: 

R, = DT, + safety level. (51) 


Finally, it should be noted that UICP assumes DLR 
items are never Mark 0 items. Therefore, only the Normal and 
Negative Binomial distributions are used to compute these 


reorder points. 
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III. THE AGGREGATE REPAIRABLES MODELS 


This chapter first discusses the proposed "Aggregate" 
Depot Level Repairables Model developed by Professor McMasters 
at the Naval Postgraduate School. Then it describes a 
simulation model which attempts to recreate the "real-world" 
repairables system. Finally, it will describe the steps of 
this model’s simulation analyses. The results of the analyses 


will be presented in Chapter IV. 


A. THE AGGREGATE MODEL 

The basic approach in the model’s development is to first 
divide the RFI inventory, inventory position, demand and 
reorder point into two separate categories, as is done in the 
current UICP integrated repairables model. The first category 
contains the inventory when the anticipated demand is for a 
unit with a returned carcass which will eventuaily be 
successfully repaired. The other category will contain the 
inventory when the carcass will be either lost or determined 
to be uneconomical to repair. 

1. Definition of terms 

Define Q, as the procurement quantity of new material 

each time a procurement order is generated, and Q, as the 


quantity inducted into a depot every time a repair order is 
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generated. Let D, be the quarterly demand rate of those units 
for which a carcass is turned in and will eventually be 
repaired, and call it the "repairable" demand rate. Let D, be 
the quarterly demand rate of those units for which a carcass 
is either not turned in or, if turned in, cannot be repaired. 
D, is called the “attrition" demand rate. Then D, and D, can 


be computed as 
D, = (CRR) (RSR) (D) = (2) (52) 


and 





” 


D, = [1 - (CRR) (RSR)]D = oo Ip. (53) 


Adding equations (52) and (53) results in 
D = dD, + D3, (54) 


where D is the total forecasted quarterly demand. Assume that 
both D, and D, are Poisson distributed. It is then also true 
that D is Poisson distributed. [Ref. 8] 

Next, let the RFI inventory position for each category 
be denoted as IP, and IP;, respectively, and let IP denote the 


aggregate RFI inventory position. It follows that 
IP = IP, + IP,. (55) 


Finally, let the maximum RFI inventory position of each 


category be denoted as SW, and SW,, respectively, and let SW 
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denote the maximum aggregate inventory position. Then, from 


equation (55), 
SW = SW, + SW,. (56) 


2. Inventory Position Ranges and Probability 
Distributions 
As demands of each type occur, the corresponding 
inventory positions decrease one unit at a time until either 
there are Q, carcasses or Q, attritions. Assume that a repair 
order or an order for procurement will be generated instan- 
taneously whenever needed. Each will then instantaneously 
increase the corresponding inventory position by Q, or Q,. 
Figure 1 shows the state transition diagram for the 
inventory position of the attrition demand category. Note 
that IP, ranges from SW, down to SW, - (Q, - 1). Because of 


the Poisson process, it follows from Reference 3 that 


Pp(SW, - w) = for w=0,1,2,...0,-1, (57) 


a 
Op 
where w represents the number of units below the maximum 
inventory position, SW,. 

Similarly, Figure 2 shows the state transition diagram 
for the inventory position of the repair demand category. 
From the same arguments used in the attrition demand case, 


Pp(SW, - v) = > for v= 0,1,2,...0,-1, (58) 
R 
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Figure 1. State Transition Diagram for Inventory Position 
of Attrition Demand. 


Inducted 





Figure 2. State Transition Diagram for Inventory Position 
of the Repair Demand. 


where v is the number of units below the maximum inventory 
position, SW,. 

In reality, the attritions from the inducted batch are 
actually not known until after the batch is inducted and 
repair is attempted. However, in the steady state situation, 
the previously inducted batch can be considered as generating 


some of the current attrition demands while carcasses are 
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accumulating for the next batch. In addition, when inducted 






carcasses are determined to be not repairable, the failure is 


usually quickly reported to the ICP’s so that attrition demand 






does not wait until the entire batch is processed through a 










depot. Thus, while D, was assumed to represent continuous 





attrition demand, in reality only that part corresponding to 






demands without carcass turn-ins is continuous. The attrition 






from repair is only approximately continuous. 






+» The Aggregate Inventory Position Probability 









Distribution 


Define x as the aggregate number of units below the 






maximum total inventory. position, SW. Then, 






SW- x = SW, - w+ SW, - v, (59) 







X=wrey. (60) 






From equations (57) and (58) the ranges of w and v are 






respectively, 








Adding (61) and (62) yields 


Osw+vs Op+ Op - 2: (63) 

or 
Os X S$ 0Op+O,- 2. (64) 
From the definition of x and equations (56), (57), 


(58), and (60), the SW term will no longer be used in the 
arguments. Rather, only w, v, and x will be considered. 
Counting from SW downward, the probability distribution for 
the total inventory position can be derived as the convolution 


of the two Uniform distributions, p, and pg. That is, 


p(x) = DY pplw)pglv = x - w), (65) 


weo0 


for x = 0,1,2,...(Qp + Q, - 2). Applying equations (57) and 
(58) to (65), subject to the inequalities (61), (62), and (64) 


results in: 





x+i1 
for 0s x<x 
O,0z : 
Min (Op, Og) for xX, CxS Xz 
p(x) = O7Qp , (66) 
Xmax + 1 - x 
at —— for x, i XS 
O;0r 2 ‘Max 
0 otherwise. 
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where x,., X,, and x, are defined as 


Xyax = Op + On - 23 
xX, = Min (Op, Og) - 1: (67) 
X2 = Xyay ~ Xy- 

The next step is to model the aggregate net inventory 
level. Unfortunately, this is analytically difficult at 
present. Arguments such as those presented on pages 28 
through 31 have been attempted. However, the assumption was 
that a batch of carcasses inducted into repair were all 
finished in a total time T,. In reality, this is not true. 
T, is the time an individual carcass is in repair, while 
carcasses are inducted serially. Thus, reality is a single- 
server queue for each batch. Therefore, this simulation model 
has been developed to provide an understanding of the nature 
of the net inventory as a function of various model parameter 


values. 


B. THE SIMULATED REPAIRABLES MODEL 

This simulation has tried to build a system which reflects 
reality as closely as possible. It takes a different view of 
the Depot Level Repairables (DLR) inventory management system 
than does UICP. It recognizes only one "aggregate" inventory 
position rather the dual view of the procurement and repair 
inventory positions of UICP. It also views attrition 


differently. Attritions occur in two different time frames. 
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First, they can occur if no turn-in accompanies the demand or 
the turn-in is lost in shipment to a supply center or other 
holding point. When this happens, it is assumed to occur 
instantly at the time of the demand. The other way to incur 
an attrition is for the unit to be condemned at the depot 
after it has been inducted for repair. Thus, the process 
generating each of these types of attrition are quite 
different. Because this simulation model can monitor 
attritions from the repair process, it allows the second type 
of attrition to occur at a different time. Once attritions 
reach a fixed level, called Q, in the model, an order is 
generated to buy the Q, quantity. Although UICP does not 
track attritions, it implicitly assumes that both types of 
attrition occur instantaneously upon demand of a unit. As 
described in Chapter II, it also uses a procurement reorder 
point to determine when to procure the next order. 

This model also takes a different view of the carcass 
repair process. UICP makes the assumption that at a certain 
point in time a batch of Q, units are inducted into the depot, 
and that all of them are repaired and returned to the supply 
system a Repair Turn-Around Time (RTAT or T,) later. This, in 
fact, does not actually happen. Carcasses are sent to the 
depot in batches. When a batch reaches the depot, it enters 
a serial queve. The first unit is examined to determine if it 
is repairable or not. If it is not, it is condemned 


immediately, and the next carcass is examined. If it is 
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determined to be repairable, it is inducted into the repair 


process and is returned as an RFI unit to the supply center a 
time T, later. This carcass passes through the first work 
station in a processing time, REP. When processing is 
complete the next carcass is examined. 

To illustrate this process, assume five units are sent to 
a depot at time T. Also, assume the third and fourth units 
are not repairable. The first unit is examined at time T and 
returned in RFI condition at time T+ T,. The second carcass 
is examined at time T + REP and returned in RFI condition at 
time T + REP + T,. The third unit is examined at time 
T + 2(REP), and condemned immediately. Because the third unit 
had been condemned, the fourth carcass is also examined at 
time T + 2(REP). It, too, is condemned immediately. Finally, 
the fifth unit is examined at time T + 2(REP) and found to be 
repairable. It is then returned in RFI condition at time T + 
2(REP) + T,. 

1. Model Assumptions 

The Simulated Repairables Mode] makes some of the same 

basic assumptions as the UICP DLR model. The model 


assumptions are listed below. 


1. Demand follows a Poisson distribution. 


2. Carcass return rate, CRR, repair survival rate, RSR, 
procurement lead time, L, repair turn-around time, T,, 
and the time it takes a repairable unit to pass through 
the first work station, REP, are known and constant. 
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10. 


11. 


12. 


Cc. 


. Demands occur at a rate of one unit per requisition. 


Carcass turn-ins occur instantaneously when an associated 
demand occurs. 


Shipments of carcasses from a supply center or other 
holding point to the depot occurs instantaneously. 


Inventory position changes with each demand, attrition, 
reorder of new material, and the induction of carcasses 
into repair at the depot. 


Carcasses are inducted in batches. Once in the repair 
system, it is assumed that the first carcass in the batch 
is immediately reviewed to determine if it can be 
successfully repaired. If it cannot, it is condemned 
instantly and the next carcass is reviewed. If that 
carcass is determined to repairable, it begins the repair 
process. The next carcass is then reviewed a REP time 
later. 


Once a carcass begins repair, it is returned to ready- 
for-issue inventory (INV) a time T, later. 


Purchases are received sequentially in the order they 
were placed. 


A purchase order is generated when the number of 
attritions, ATTR, reaches a specified order quantity, Q,. 


A repair order is generated when the number of carcasses 
accumulated in a repair queue, RQUE, reaches a specified 
value Q,. 


The safety level is defined as the average net inventory 
when a procurement of new material is received or a 
repaired item is returned to the supply system. 


DISCRETE EVENT MODELING 


The simulation software package used for this simulation 


model is called SIGMA, the SImulation Graphical Modeling and 


Analysis system, developed by L. Schruben, Cornell University, 


and distributed by Scientific Press, South San Francisco, 
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California. It was originally developed to study the dynamic 


behavior of systems; i.e., how systems change over time, and 
is ideally suited for simulating the Navy’s repairables 
inventory management process. It is an interactive graphics 
approach for building, testing, and experimenting with 
discrete event simulations on personal computers. A powerful 
feature of SIGMA is that it can automatically translate a 
simulation model into portable source codes for Pascal and C, 
and can be compiled and run on a wide variety of computers. 
(Ref. 9] 
1. Discrete Event Simulation Modeling Terms 

To discuss how discrete event simulation works, 
certain terms need to be described. A "system" is defined as, 
",..a collection of entities that interact with a common 
purpose according to sets of laws and policies." [Ref 9] 
Systems are characterized by their function or purpose such as 
a transportation system or an inventory control system. The 
"entities" are the elements which are used to perpetuate or 
"fuel" the system. Examples of entities are components on an 
assembly line system or variables in a system of equations. 
They are described by their “attributes" or characteristics. 

A discrete event system is one in which changes in the 
system occur at particular instances of time. The "state" of 
a system is simply a description of that system. In a 


computer simulation of a discrete event system, the state is 
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defined by the numerical values of its variables (referred to 
as the state variables) and the schedule of future “events” 
(state changes within the system). Finally, “laws and 
policies" determine the conditions under which future events 
will occur, the time they are scheduled to occur, and the 
priority of execution between two or more events scheduled at 
the same time. 
2. Discrete Event Simulation 

In a discrete event system simulation time is advanced 
in discrete steps until it reaches a value at which the state 
is scheduled to change. At this point, one or more variables 
change value, and one or more events may be scheduled or 
canceled. An excellent example is the Navy’s repairable 
Management system. The system is changed when a demand of one 
unit occurs. This event causes the inventory position and the 
on-hand inventory to reduce their respective values by one 
unit. Other events may be scheduled, such as another demand, 
a reorder, or an event which determines if a failed carcass is 
turned in to the supply system. These events can only occur 
if certain conditions of the new state are present. They can 
be scheduled (or canceled) after a fixed, or random, time 
period interval following the current time. 

In SIGMA a main controller executes event routines. 
This controller operates from a master appointment list of 


scheduled events referred to as the "future events list" 
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because, at any given time, it contains all events which have 


been scheduled and the times they are to be executed. During 
a simulation, the main control program will advance the 
simulated time to the time of the next scheduled event, remove 
the event from the list and execute it. The main control 
program continues calling the next scheduled event until a 
condition for stopping the simulation is met. SIGMA will stop 
the simulation either after the simulation clock reaches a 
user designated time, or after a user designated event occurs 
a given number of times. This parameter is established in 
SIGMA’s RUN menu just prior to executing the simulation. A 
general flow chart summarizing this process is illustrated in 
Figure 3. 
3. Modeling with Event Graphs 

The major thrust of SIGMA is to model or simulate a 
discrete event system using an event graph. As mentioned 
above, the three elements of a discrete event system are the 
state variables, the events that change state variables, and 
the relationships between events. An event graph organizes 
these into the simulation model. In SIGMA, events are 
represented by vertices, and the relationships between the 
events are represented by directed edges. A picture of the 
event graph used in this thesis research is shown in Figure 4. 

To develop a simulation model, the three elements need 


to be defined in mathematical terms. SIGMA provides the means 
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Pigure 3. Main Event-Scheduling Algorithm. 


to define the terms using "pop-up" menus. State variables can 
be single variables, arrays, or matrices. They are defined as 
either real or integer. The events (vertices) contain 
formulas which determine how state variables change their 
values. The edges represent "laws and policies"; i.e., the 
conditions, time delay, and priority, under which one event 


schedules (or cancels) the connecting event. These conditions 
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Figure 4. Event graph of THSVER4.MOD 


are implicit "If.., then" statements which either schedule the 
connecting event if the statement is true or do nothing if the 
statement is false. 

In SIGMA, the default edge condition is "“l==#1". (The 
use of the double equal sign is peculiar to SIGMA and is only 
used in condition statements on the edges.) The edges 
meeting this condition are labeled as "unconditional" edges 
because they automatically schedule the connecting event. 
The time delay can be a number, a variable, a random number, 
or a function of any of these. SIGMA has its own random 
number generator which can be called from either an event or 


an edge. It generates an uniform random variable, RND, whose 
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range is the set of all real numbers between zero and one, 
i.e., 0.0 < RND «< 1.0. 

The use of priorities on an edge is the way in which 
SIGMA executes events scheduled to occur at the same time. 
When the scheduling condition is met by two or more out-edges 
which have the same time delay, SIGMA executes the connecting 
event with the highest priority first. The priority range is 
from one to ten with one the highest and ten the lowest 
priority. They are assigned to each event by the user. 

Descriptions of the state variables, events, and edges 
for the three models used in this thesis research 
(THSVER4.MOD, THSVER7.MOD, and THSVER8.MOD) are provided in 
Appendices A, D, and F, respectively. A short synopsis of the 
relationships of each vertex with its out-edges and connecting 
vertices is automatically generated by SIGMA. These synopses 
are provided in Appendices B, E, and G, respectively, for each 
of the three models. A detailed description of THSVER4.MOD, 
which illustrates these relationships, is discussed in section 
D below. 

4. Output 

Obtaining data from a SIGMA simulation is not a 
difficult task. The user simply designates the state 
variables and events which need to be "traced". Once control 


passes to a traced event and state changes occur, SIGMA will 


write the name of the event, the number of times the event 











occurred, and the value of all traced variables to an external 


file. An example of an output file is provided in Appendix C. 


D. THE MODEL 
1. Model Description 
The model begins with the RUN event. Its purpose is 
to initialize various state variables. Variables not given 
specific values by the user in this event are automatically 
initialized as zero. [Ref. 9:p. 30] For the repairables 
simulation model analyzed in this thesis, certain state 


variables are never changed. Their names and values are: 


L = 8.2 quarters 
T2 = 1.3 quarters 
D=9 units per quarter 
OP = 10 units 
OR = 5 units 
TIME = 850 quarters 


SW =72 units. 


Three other variables initialized in the RUN event will be 
changed from one simulation run to another. These variables 
are CRR and RSR, which range from 0.0 to 1.0 in increments of 
0.2, and REP, which ranges from 0.0 to the value of T2 (1.3) 
in increments of 0.25 quarters. 

Following initialization of the model variables, the 
RUN event schedules the INIT event. In the INIT event, IP and 
INV are set to their maximum values. From Figure 4, there are 
two edges originating from INIT connecting with DMND1 and 


STATS. As noted in Appendix D, both are "unconditional" 
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edges. The edge going to STATS has a time delay of TIME 
quarters which was used in this model to stop the simulation. 
The user must designate STATS as the event which stops the 
simulation, and must also specify that it occurs only once. 
This ensures that the simulation will stop when the time of 
the simulation reaches TIME. 

The event DMND1 schedules all demands. Thus, it 
unconditionally passes control to the DMAND event. Once 
control passes to DMAND, the state immediately changes by 
reducing the state variables IP and INV by one value, (i.e., 
mathematically, INV = INV - 1 and IP = IP - 1), and by 
assigning the state variable X[1] the first value of the 
uniform random variable, RND, to be used in the decision 
associated with the implicit Bernoulli trial from the 
conditional statements of two of DMAND’s out-edges. 

Once changes occur, the main control program schedules 
events Sasea on the conditions described on the edges. There 
are three out-edges from DMAND connecting with events CARTN, 
DMND1, and ATTRI, and denoted by edges (1), (6), and (12), 
respectively, in Appendix A. Edges (1) and (6) model a 
Bernoulli trial to determine if a failed carcass is returned 
to the supply system when a demand occurs. If the value of 
X{1] is less than or equal to CRR, then the carcass has been 
returned and CARTN will be scheduled at a time equal to the 
current time. If X[1] is greater than CRR, then the carcass 


was not returned and the condition of edge (6) is met. ATTRI 
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is, therefore, scheduled to occur immediately. If the carcass 
is returned, the variable, RQUE, will increase by one unit. 
If not, the variable, ATTRI, will increase by one unit. 

The other directed edge from DMAND connects with 
DMND1. Its condition is "l==1", making it an unconditional 
edge. In other words, the event DMND1 will always be 
scheduled whenever DMAND occurs. The time delay of scheduling 
DMND1 is (-1/D)*ln(RND), which is the exponentially 
distributed random time between demands and is a consequence 
of the Poisson process. Therefore, DMND1 will be scheduled at 
the current time plus the delay. In the case where the time 
delay is equal to zero, -there will be two events, DMND1 and 
either CARTN or ATTRI, scheduled at the same time. However, 
the edges scheduling CARTN and ATTRI have been assigned a 
priority of three by the author, and the edge scheduling DMND1 
has been assigned a priority of five so that CARTN or ATTRI 
will be executed before DMND1. 

Demands continue to occur in the above manner until 
the simulation stops. As demands occur, the variables ATTR 
and RQUE increase in size. When either reaches a given level, 
the respective procurement and repair processes begin. 

2. Repair Process 

Once RQUE reaches the value QR, the batch of carcasses 

in RQUE are inducted into the repair process. As a 


consequence, IP and QUEREP are increased by QR, and the event 
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RPORD is immediately scheduled. This event is the focus of 
the repair process because it is the event which determines 
whether a carcass can be repaired or must be condemned. The 
event RPORD, coupled with the COUNT event, will actually 
analyze each carcass sequentially (but at the same simulated 
clock time) and make a decision on each carcass. Then, as a 
result of the edge time delays, each carcass’ delay in 
entering attrition or repair categories will be assigned a 
repair processing time within the range of zero to 
(QR - 1) * REP. 

Just prior to the occurrence of RPORD, the CARTN event 
will have just increased IP and QUEREP by QR, and set the 
variable COUNT equal to zero. Then the event RPORD is 
executed in a string of QR distinct occurrences at the same 
simulated clock time. As this sequence proceeds, RQUE empties 
one at a time as the string of Bernoulli trials is executed 
and the decision is made to either repair or condemn the 
carcass. If the carcass can be repaired, the event DEPOT is 
scheduled. If it cannot, the event ATTRI is scheduled. The 
time delay effect of REP is then accessed through the COUNT 
event. 

The time delay on the edges from RPORD to ATTR1 and 
DEPOT models the repair processing time. As each carcass is 
analyzed, a time delay is assigned based on the product of REP 
and the variable COUNT. After the first carcass is examined 


in RPORD, the time delay in scheduling either ATTRI or DEPOT 
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is COUNT * REP = 0, since initially the variable COUNT equals 
zero. The return to RFI condition for this carcass, then, 
would be the current time plus T2. At the same time, but 
after either ATTRI or DEPOT is scheduled, the event, COUNT, is 
scheduled and executed. The only state change in this event 
is to change the variable, COUNT. If the first carcass is 
repairable, then the variable COUNT increases by one. If not, 
COUNT will remain zero. If there are any more carcasses in 
the batch, (i.e., if RQUE is greater than zero) then control 
passes immediately back to the event RPORD where the next 
carcass is processed. (Notice, no simulated clock time has 
expired.) 

Since only one carcass has been processed at this 
point (and assuming QR > 0), RPORD examines the second 
carcass. RQUE is reduced by one and another random number is 
generated for the Bernoulli trial on the edges leading to the 
events ATTRI and DEPOT. The time delay will now have been 
affected by the first carcass. If the first carcass was 
repairable, then COUNT would have been increased to one, and 
the time delay would be REP times the quantity COUNT, or a 
total of REP quarters. The carcass would then be returned to 
RFI condition at the current time plus REP plus T2. If the 
first carcass was condemned, COUNT would not increase, and the 
time delay for this carcass would be zero. This process 
continues until RQUE equals zero and control passes to the 


next scheduled event. 
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An illustration of how this simulation models the 
repair process is presented here using typical values for the 
variables. Assume RSR equals 0.8, REP equals 0.25, QUEREP 
equals zero, IP equals 60, and the simulated clock time is 18. 
Also, assume the constant variables mentioned above in Section 
D.1 retain their designated values. 

Suppose at time 18, the variable RQUE reaches QR = 5 
in event CARTN. This will cause the variables QUEREP and IP 
to be increased to five and 65, respectively. The condition 
on the out-edge to RPORD is met; therefore, RPORD is 
immediately scheduled and executed. This causes RQUE to be 
reduced by one unit, to four. Suppose that X[0)’s random 
variable value is found to equal 0.6. Since X{0]) is less than 
RSR = 0.8, the edge condition from DMAND to DEPOT is met and 
event DEPOT is scheduled to be executed at time 18 + 0 = 18 
since COUNT is zero and the product of COUNT and REP is zero. 
Control moves to the event, DEPOT, since its edge has a lower 
numbered priority than the edge leading to the event COUNT. 
There are no state variable changes here, but the event, RPAIR 
is scheduled to occur in T2 time units (i.e., the edge from 
DEPOT to RPAIR simulates the repair turn-around time). 
Control will pass to RPAIR at the current time plus T2 and the 
first carcass will be simulated to return to RFI condition at 
time 18 + 0 + 1.3 = 19.3 quarters. 

Immediately following this execution of DEPOT, the 


event COUNT is executed. Since X[0] is less than RSR, the 
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variable COUNT increases to a value of one and schedules RPORD 
with no delay. (At this point in the example, control in the 
repair process will fluctuate between events RPORD and COUNT 
until RQUE becomes zero.) RQUE is reduced by one more unit to 
three and another value for X[{0]) is generated, assumed here to 
be 0.4. Because X[0] is less than RSR, DEPOT will be 
scheduled with a time delay of REP * COUNT = .25 * 1 = 0.25 
quarters. Once DEPOT is executed, the event RPAIR will be 
scheduled at a time delay of T2. This second carcass will be 
returned to RFI condition at time 18 + .25 + 1.3 = 19.55 
quarters. 

The event, COUNT, will be executed again, increasing 
the variable COUNT by one unit to a total of two, and for the 
third time in this sequence, RPORD is executed. RQUE is 
reduced by one unit to two, and a new value for xX[0] is 
generated, assumed here to be 0.85. Since X[0] is now greater 
than RSR, the event ATTR1 will be scheduled to occur at the 
current time plus REP * COUNT = 0.25 * 2 = 0.5 quarters, or at 
clock time 18 + 0.5 = 18.5 quarters. At that time, the repair 
attrition will be simulated, reducing IP and QUEREP by one 
unit. 

The event COUNT is again executed at clock time equal 
to 18. However, since X[0] was greater than RSR, COUNT will 
not be increased, retaining its value at two. RPORD will then 


be executed. 
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This is the fourth time in this sequence that RPORD 
has been executed. RQUE is now reduced one unit to the value 
one and another value for X[0) is generated and assumed to be 
0.92. As in the previous case, ATTR1 is again scheduled to 
occur at time 18.5 quarters, resulting in the reduction of 
QUEREP and IP by one unit. 

The event, COUNT, is executed a fourth time at time 18 
quarters, and again, because this carcass was not repairable, 
the variable COUNT remains at two, and RPORD is executed. In 
RPORD, RQUE is reduced by one unit to zero and the next value 
for X[{0]) is generated; say, 0.2. This implies that DEPOT will 
be scheduled at time 18 + .5 = 18.5 quarters and it will 
schedule the event RPAIR to occur at time 18.5 + 1.3 = 19.8 
quarters. 

The event COUNT is scheduled and executed once again 
(still at time 18 quarters), increasing the variable COUNT to 
three. This event becomes meaningless, however, as RPORD will 
not be scheduled since RQUE is equal to zero and no longer 
meets the condition on the edge from COUNT to RPORD. Thus, 
all that is left in the repair process for this batch is the 
occurrence of the scheduled repairs and attritions. 

To complete the repair process, every time RPAIR 
occurs QUEREP is reduced by one unit and the event SHIP1 is 
immediately executed, allowing the net inventory, INV, to be 
increased by one unit. Likewise, every time ATTR1 occurs, 


QUEREP decreases by one unit and the event ATTR2 is executed 
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increasing the attrition variable, ATTR, by one unit. As 
explained previously in this example, two attritions occur 
simultaneously at time 18.5 quarters, and three carcasses will 
be returned to the supply system in RFI condition at times 
19.3, 19.55, and 19.8 quarters. 

One of the purposes of THSVER4.MOD was to record the 
inventory position, IP, each time it changed its value. IP 
changes only in events DMAND, ATTR2, CARTN, and ORDER. 
However, in simulations where CRR is greater than zero and QR 
is greater than one, it does not always change in the event 
CARTN. IP changes in the event CARTN only when the value RQUE 
increases to the value QR. This leads to difficulties when 
recording only IP changes. Tracing CARTN would be impractical 
because each time CARTN is executed, IP would be recorded 
whether it had changed or not. 

To avoid this problem, a "dummy" event, named INPOS, 
was created. INPOS was connected to the event COUNT where its 
condition was "RQUE == 0" and its delay time was also zero. 
This resulted in control passing to the event INPOS only after 
all carcasses had been examined. Thus, tracing INPOS for IP 
would result in the recording of IP only when the change in IP 
is due to the induction of QR carcasses. 

3. Procurement Process 
The procurement process is much less complex than the 


repair process. procurement is generated whenever the 
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inventory system loses a given number of assets from 
attrition. Attritions occur whenever an asset fails to enter 
the repair process (e.g., lost in shipment or surveyed by the 
end user) or is condemned at the repair facility. Attrition 
due to loss or survey is modeled as a Bernoulli trial in event 
DMAND, based on the carcass return rate, CRR. Attrition due 
to condemnation is modeled as a Bernoulli trial in the event 
RPORD, based on the repair survival rate, RSR. Every time a 
carcass is lost or condemned, the variable ATTR is increased 
by one unit. Once ATTR reaches the value QP, an order is 
generated, IP is increased, reflecting the order, the variable 
QUEORD is also increased by QP, and ATTR is reset to zero. 
The receipt of the order is scheduled for L = 8.3 quarters 
later in time. INV will increase by QP at that time, while 
QUEORD will decrease by QP. 
4. Steady State Determination and the Termination Process 
As in many time series simulations, this model needs 
to run for a period of time to "settle" into a steady state. 
Appendix H contains various graphs of the net inventory versus 
time with widely different values of CRR and RSR. From a 
visual inspection it appears that steady state is achieved 
after approximately 10 quarters. As a precautionary measure 
in this study, an additional ten quarters were allowed to pass 


before data were recorded for the analyses to be presented in 
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the next chapter (i.e., the first 20 quarters of data were 


ignored) . 

To determine when to stop the simulation runs, the 
limits of the APL2 editor, called Editor 3, needed to be taken 
into account. This editor is provided by International 


Business Machines as part of the PC version of its APL 232 


software. It was used initially to edit each simulation 
output. Originally, a simulation was stopped after 360 
simulated quarters. However, when graphing the resulting 


simulated data into a histogram, the shape of the graph 
appeared "ragged". To "smooth" the shape, the time of the run 
was extended to 1,200 quarters, however, Editor 3 was unable 
to capture all the data when the edited version was saved. It 
was discovered that if the simulation stopped after 850 
simulated quarters, the edited data would be retained in full. 
Therefore, 850 quarters were used as the ending time of each 
simulation. (Later, the MS-DOS editor, Edlin, was suggested 
as an alternative editor for the purpose of this research. It 
was found to be a better tool than any other editor available. 
It can easily edit the 1,200 quarter data sets. In fact, it 
was used to edit the largest single file (39 megs of data) 
generated from a batch run. Unfortunately, by the time this 
discovery was made, the vast majority of the simulation runs 
had been made with the output edited using the time of 850 
quarters. Thus, to remain consistent, 850 quarters was used 


as the ending time for all simulation runs.) 
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5. Revisions to the Model 

SIGMA, as noted, is a powerful simulation software 
package. It has two major features which make it powerful. 
First, it has the ability to read data from an external file. 
This feature is important when running batch simulations 
because variables can easily be assigned a range of values in 
an external file. Second, it has the ability to automatically 
reset all variables to zero, reset the random number seed to 
the original value, cancel all events on the future events 
list, and restart the simulated clock time by using the SIGMA 
function, "SET { }", on an edge. The combination of both of 
these features allows the running of batched simulations. 

Initially, individual simulation runs of THSVER4.MOD 
were needed to graph and examine the resulting histograms 
associated with inventory position and net inventory for each 
combination of the ranges of CRR, RSR, and REP listed above. 
These histograms were used to identify the nature of the 
probability distributions of the Inventory Position and Net 
Inventory and to examine the effect of the three variables on 
them. Later analyses did not require these graphs, and a 
model, THSVER7.MOD, designed for batch runs was developed. 
This model includes an exact duplicate of THSVER4.MOD and, in 
addition, the SET { } function is introduced on an out-edge 
from event STATS to a new event, RINIT. When the event STATS 
occurs, the introduction of SET { } on the edge automatically 


resets and restarts the model as noted above. RINIT 
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reinitializes the variables L, T2, QP, QR, and TIME to their 


original values. The out-edge from RINIT to INIT is 
unconditional with no delay. Thus, INIT is scheduled and 
reads the next set of values of D, CRR, RSR, and REP from an 
external file, THS.DAT. (For each simulation run of this 
thesis, the value, D, remained at nine units per quarter. 
Future thesis efforts can use this feature to study the 
effects of changing D.) 

This model was first used to determine inventory 
safety levels for all values of CRR, RSR and REP. As 
mentioned in the model assumptions, the safety level is 
defined as the average net inventory at the time either a 
repaired carcass returns to the supply system in RFI 
condition, or a procurement order is received. The model also 
introduced two new variables, SAFTY and COUNT[{1). Both were 
used in the events where INV changed as a result of these 
types of receipts (i.e., in the events SHIP1 and RCVD2). In 
each event, SAFTY summed the value of net inventory; i.e., 
SAFTY = SAFTY + INV, while COUNT[1] accumulated the total 
number of times net inventory increased. When 850 quarters 
was reached, the event STATS computed the average safety level 
using the formula, SAFTY = SAFTY/COUNT[1]. Appendices D and 
E provide the details of the THSVER7.MOD model. 

To study the average time-weighted values of net 
inventory, on hand inventory, and backorders, an additional 


model was developed. Using the batch model discussed above, 
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the new model, THSVER8.MOD, introduced an array of variables 
to accumulate this data. A six-variable array, named Cl[j], 
was required as SIGMA limits the number of variable names to 
20 and assigning six additional names rather that using an 
array would have resulted in 24 names, exceeding the 


limitation. The elements of C[j] are described as follows: 


C[1}: the time of the last net inventory change; 


4 


2. C(2): the time of the current net inventory change; 

3. C{3): the difference between C[2]) and C[1]. It is the 
time which elapsed before the net inventory changed 
value; 


4. C[4]: accumulates the total time-weighted segments of net 
inventory at any time INV changed value for the run; 


5. C(5]): accumulates the total time-weighted segments of INV 
whenever INV was greater than or equal to zero at any 
time it changed value during the run; 

6. C({6]: accumulates the total time-weighted segments of INV 
whenever INV was less than zero at any time it changed 
value during the run. 

Because INV changes values in only three events 
(DMAND, SHIP1, and RCVD2), extracting the time-weighted data 
was not difficult. The relationships between the above 
variables were formulated in the event just prior to the 


change in INV (i.e., in any of the events DMND1, RPAIR, and 


RCVD2). The relationships are as follows: 


1. C1) 


C{2); sets C[1] to previous clock time. 


2. C2) 


CLK; sets Clz) to current clock time. 








3. C{3]) = (C{2]) - C[1]), if CLK 2 20; computes the time 
weighting factor for time intervals after the 20th 
quarter. 


4. C{4) = cC{4) + C[3]) * INV; computes the total time 
weighted value of INV. 


5. c{5]) = C(5) + C[3] * INV, if INV 2 0; computes the total 
time-weighted value of O/H inventory. 


6. cC{6) = C(6) + C{3) * INV, if INV < 0; computes the total 
time-weighted value of backorders. 


In the STATS event, the averages were written using the 


following formulas: 


c{4] = c[4]/830, 
c{5] = c{5)/830, 
= C[6])/(-830). 


c[6] 


Appendices F and G together provide a detailed description of 


THSVER8 .MOD. 
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IV. RESULTS 

The purpose of this chapter is to discuss the significant 
results of the simulation study. The first section provides 
a comparison of the simulated results with the theoretical 
trapezoidal distribution of the time-weighted inventory 
position presented in Chapter III. The next section discusses 
the safety level analysis and compares the safety level 
calculated by the UICP model with the proposed method of 
defining the safety level described in Chapter III. Finally, 
the third section presents a discussion of the results of the 


simulated aggregate net inventory distribution. 


A. TIME-WEIGHTED INVENTORY POSITION PROBABILITY DISTRIBUTION 

To validate the probability distribution of the time- 
weighted inventory position which was described in Chapter 
III, a set of batched simulations were run using THSVER8.MOD. 
The variable IP was traced on all the events where its value 
changed; i.e., DMAND, ATTR2, INPOS, and ORDER. There were a 
total of 181 runs, one run with CRR equalling zero and the 
remaining 180 were the combinations of varying CRR from 0.2 to 
1.0 in steps of 0.2, RSR from 0.0 to 1.0 in steps of 0.2, and 
REP from 0.0 to 1.3 in steps of 0.25. Then, for each set of 
parameters, the theoretical trapezoidal probability 


distribution discussed in Chapter III was compared to the 
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results. A Chi-square goodness of fit test was attempted to 
determine if the simulation results would validate the 
theoretical distribution. For this test the "observed 
frequencies" resulting from the simulation model were 
generated by determining the total length of time that each 
value of the inventory position occurred during the 830 
periods. Thus, the frequency values are not integer. The 
Chi-square test procedure assumes integer frequency values. 
As a consequence, this use of the Chi-square goodness of fit 
test "stretches" the theory and is not exact. However, no 
other goodness of fit test exists for this type of frequency 
determination. 

The results of the Chi-square test for "goodness of fit" 
strongly supports the theoretical model. Only the simulation 
runs where CRR equalled zero, or when CRR equalled one and RSR 
equalled zero or one, did the results fail the test (the Chi- 
square statistics were greater than 500 resulting in a p-value 
of zero). For these cases, the inventory system reduces to 
the consumable model of Reference 3, and as such, the 
inventory position is Uniformly distributed. 

The remaining simulation runs resulted in a Chi-square 
statistic for each case well below the value of 27.7 or 22.4 
which are the Chi-square critical values at 0.01 and 0.05 
significance levels, respectively, with 13 degrees of freedom. 
(The 13 degrees of freedom were derived from the fact that 


there were 14 possible values IP could take and none of the 
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expected frequencies were computed using any estimated data. 


Therefore, the degrees of freedom were 14 - 1 = 13.) This 
leads to the non-rejection of the null hypothesis that the 
empirical distribution of the simulated data follows this 
trapezoidal distribution. The test results are presented in 
Appendix I. 

An important result of this analysis is that while the 
Chi-square statistic changes with REP, it is not statistically 
Significant. As described in Chapter III, the theoretical 
distribution assumes REP equals zero in deriving the formulas 
for D, and D,. The variable REP introduces a delay in 
attrition accumulation. However, for these simulations, REP 
does not seem to have an effect on the inventory position 
distribution. Even when REP has its strongest influence on 
the inventory system (i.e., when all the carcasses are 
repaired (RSR = 1.0) and REP equals the repair turnaround-time 
(T,) of 1.3 quarters) the Chi-square statistic did not change 
significantly. This would suggest, then, that the model 
discussed in Chapter III is robust. However, further analysis 
involving changes in Q,, Q,, D, lL, and T, are needed before 


this claim can be fully justified. 


B. SAFETY LEVELS 
As noted in assumption (12) of Chapter III, safety level 
is defined as the expected value of the net inventory at the 


instant a procurement order or successfully repaired carcass 
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is received at a stock point. To find this aggregate safety 


level for the inventory system being modeled, THSVER7.MOD was 
used for a set of batched runs which varied the parameters 
CRR, RSR, and REP as mentioned above. The resulting simulated 
safety levels are found in Appendix J. 

A "theoretical" integrated repairables model safety level 
using the Navy’s UICP formulas was also computed using the 
Same range of values for CRR and RSR. These values take into 
account the dual inventory positions established in UICP as 
described in Chapter III, section A. As discussed in Chapter 
III, equation (56) shows the maximum inventory position is, in 
fact, equal to the sum of the maximum inventory positions of 


the repair and procurement processes: 
SW = SW, + SWp. (56) 


Under the UICP model, the maximum values of SW, and SW, occur 
at the instant a repair order or a procurement order is 
generated, respectively. The UICP model describes the maximum 


inventory positions for each part as 


SW, 


2, + SAFETY LEVEL + Q, (68) 
and 


SW, 


Z + SAFETY LEVEL + Qp, (69) 


where the safety level is the same in both cases for the 
integrated repairables model, as explained in Chapter III 


(Ref. 6]. 
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By substituting equations (68) and (69) into equation (56) 


one can compute the UICP total safety level, which is the sum 


of the two safety levels shown above, as: 
TOTAL SAFETY LEVEL = SW- (Z + Z, + Op + Op): (70) 


when SW, Q,, and Q, are fixed, and Z and Z, are computed using 
the formulas from Chapter II. These results are compared with 
the simulated results for REP = 0.0 in Table 1. The REP = 0.0 
case results in all repaired carcasses in a batch being 
returned after a time interval of T,. Interestingly, in 
almost all cases, the simulated safety levels exceeded the 
UICP safety levels. They even exceeded the values of the UICP 
safety levels divided by two when RSR 2 0.2. This suggests 
that the UICP does not accurately reflect the "real world" of 
the depot level repairables inventory system and, in fact, 


would understate the actual aggregate safety levels. 


C. TIME-WEIGHTED NET INVENTORY 
1. Net Inventory Distribution 
The Hadley-Whitin Consumable item inventory system 
described in Chapter II and Reference 3 allows backorders and, 
in particular, incorporates the expected number of unit years 
of shortage (i.e., time-weighted backorders) incurred per year 
in the backorders cost term of the Total Variable Cost 


equation. Of course, the expected time-weighted on-hand 
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TABLE 1. COMPARISON OF UICP COMPUTED SAFETY 
LEVELS WITH SIMULATED SAFETY LEVELS. 


inventory is included in the holding costs term. These terms 
affect both the determination of the Order quantity and 
Reorder Level. 

The expected unit years of on-hand and shortages can 


be obtained from the probability distribution of time-weighted 
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net inventory (i.e., w(x) for all -w < x < x). The formulas 
for the time-weighted expected on-hand inventory (O/H) and 
time-weighted expected number of backorders (B/O) are given by 
equations (71) and (72), respectively: 


Sw 
o/H = So xv (x) (71) 


x20 


° 
B/O = -)0 xV(x) . (72) 


For the consumable model, when demand is generated by a 
Poisson process, lead time is deterministic, and the lead time 
demand is large, the probability distribution of demand during 
lead time can be approximated using the Normal distribution. 
(Ref. 3:p. 192] The probability density function for the net 
inventory is then, 
meee ra aera a (eee cians |: (73) 
for all -~ s x < SW, 

where Q is the order quantity, r is the reorder level, yp is 
the mean lead time demand, o7 = yu, and (x) is the 
complementary cumulative distribution for the Standard Normal 
function. Thus, ¥(x) is not Normally distributed. 

Equation (73) can then be used to compute the expected 
on- hand (E(0/H)) and expected backorder (E(B/O)) quantities. 
Hadley and Whitin derived the equation for expected on-hand. 


{Ref. 3:p. 194) This equation is given as 
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Xmax 


E(O/H) = [ x00) 
0 


(74) 
= 2+ r-p + £(B/0). 
Because 
E(N/I) = E(O/H) - E(B/O), (75) 
the expected net inventory can be computed as 
E(N/I) = g +r-4p (76) 


by substituting equation (74) into (76). 

Equation (76) was then used to calculate the 
theoretical mean of net inventory for the two limiting cases 
(i.e., when CRR = 0.0, or when CRR = RSR = 1.0 and REP = 0.0), 
for the purpose of comparison with the simulated results. The 
values used in equation (76) when CRR = 0.0 (the completely 
consumable case) were 

Q=0,=10 

I = Xyy~ Op= 72-10 =62 

p=D*L=9*8.2=73.2. 
Therefore, E(N/I) = -6.80. When CRR = RSR = 1.0 and REP = 0.0 
(the completely repairable case), the values used in equation 
(76) were 

O=0p,=10 


ZT = Xyay~ Op = 72-5 =67 
peD*T, =9*1.3=11.7. 


Therefore, E(N/I) = 57.8. 
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a. Simulated Expected Value of Net Inventory 


The expected values of the simulated time-weighted 
net inventory for the various values of the parameters CRR, 
RSR, and REP were generated using THSVER8.MOD as_ the 
simulation vehicle and are presented in Appendix K. The 
resulting values for the limiting cases were then compared 
with the calculated ones using equation (76). When CRR was 
zero, the value of the simulated expected or mean net 
inventory was -6.92 while the calculated theoretical value was 
~6.80; when CRR and RSR were one and REP was zero, the 
simulated mean net inventory was 58.1 while the calculated 
value was 57.8. This validates the simulation model results 
for these cases. 

b. Simulated Net Inventory Distribution 

Net inventory distributions for various values of 
CRR, RSR and REP were generated and graphed as histograms. 
The initial graphs appeared bell-shaped (see Appendix L), so 
a Normal probability distribution was "fitted" using the AGSS 
software. The Chi-square results are displayed in Table 2. 
Rather than showing the critical Chi-square values for a = 
0.05 and 0.01, the AGSS software provides the alpha value 
which would be needed to give the computed Chi-square 
statistic. An a value of 1.0 means that virtually all the 
Chi-square distribution for the given degrees of freedom was 


to the right of the computed statistic. Thus, the probability 
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TABLE 2. SUMMARY OF CHI-SQUARE STATISTICS FOR 
FITTING NORMAL CURVE TO SIMULATED TIME- 
WEIGHTED NET INVENTORY DISTRIBUTIONS. 


of making a Type I error is virtually zero in accepting the 
null hypothesis that the data is Normally distributed at any 


@ level. The lowest value of a in Table 1 was 0.20951, which 
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means that approximately 21% of the area of the Chi-square 
distribution was to the right of the computed statistic. 
Thus, the null hypothesis is again strongly accepted at the 

a = 0.05 and 0.01 level since it would be accepted for any 
significance level less than 0.20951. 

In conclusion, the results of Chi-square goodness 
of fit tests in all cases suggest that one cannot reject the 
null hypothesis that the net inventory distribution fits the 
Normal distribution. 

What does this conclusion imply about’ the 
distribution given by equation (73) for the limiting case of 
CRR = 0.0 (the completely consumable case)? As noted earlier, 
equation (73) is not a Normal distribution. It would only be 
Normal if Q = 1.0, and that was not the case in the simulation 
since Q, = 10.0. Perhaps it is close to Normal because of the 
parameter values selected for the simulation. At this time 
nothing further can be concluded. However, the result 
suggests that further analysis of that formula is appropriate 
to see under what general conditions it can be approximated by 
a Normal distribution. 

2. Sensitivity Analysis of CRR, RSR, and REP. 
a. Mean and Variance as a function of CRR, RSR, and 

REP 

Assuming, from the above results, that net 


inventory is Normally distributed, the next step is to examine 
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the values of the two distribution parameters, the mean and 
the variance, as a function of the parameters which were 
varied; namely CRR, RSR, REP. The data was generated using 
THSVERS.MOD. The values of the means and variances of the 
simulated distributions for net inventory are summarized in 
Appendix K. Using the data from Appendix K, graphs were first 
produced for the distribution means; graphs for the variances 
were then added. These graphs are paired in Appendix M. In 
these graphs, one parameter is fixed; the second parameter 
represents the independent variable for plots of curves for 
fixed values of the third. Each line or curve is identified 
by an integer which represents a unique value of this third 
parameter. When the third parameter is REP, a line or curved 
identified with a "0" means REP = 0.0; "1" means REP = 0.25; 
"2" means REP = 0.0; "3" means REP = 0.75’ "4" means REP = 
1.0; "5S" means REP = 1.3. When the third parameter is CRR or 
RSR, the integer identifying the line or curve represents the 
value of CRR or RSR where a "0" means 0; "1" means 1.0; "8" 
means 0.8; "6" means 0.6; etc. 

In addition to these graphs, because the UICP model 
uses the product of CRR and RSR to forecast attrition demand, 
plots of the mean and variance of the net inventory 
distribution for various values of the product of CRR and RSR 
were made. These graphs are provided in Appendix N. 

The first 23 graphs of Appendix M examine the 


influence of REP for fixed values of CRR and RSF. As REP 
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increases, the mean of the distribution for net inventory 
decreases because the carcasses of an inducted batch are 
delayed longer before returning to RFI, and its variance 
increases because of the same delay. No delay occurs, of 
course, when REP = 0. The graphs show a completely linear 
relationship between REP and the distribution mean and a 
fairly linear relationship for the variance as a function of 
REP. 

The REP curves for the different RSR values show 
that as RSR increases so does the mean of the net inventory 
because the repair process is providing more and more of the 
RFI units, and it nas a much shorter "lead time" (repair turn- 
around time) than procurement. The variance of net inventory 
decreases for the same reason. The REP curves for different 
CRR values show that as CRR increases for RSR 2 0.4 the net 
inventory mean and variance behave similarly for the same 
reason. 

In the next 22 graphs of Appendix M_ the 
relationship between RSR and the distribution mean and 
variance appears to be almost linear except in the case where 
CRR = 1.0 and RSR approaches 1.0. As noted above, when RSR 
increases so does the mean, while the variance decreases. 
Interestingly, the variance curves for various REP values are 
quite close. A review of the graphs for both the mean and 
variance when REP is fixed shows the crossing of RSR curves 


for various CRR values when RSR is less than 0.4 and 0.6 for 
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the mean and variance of the net inventory, respectively. 
This CRR effect is examined next. 

The final 24 graphs of Appendix M show plots of the 
distribution mean and variance as a function of CRR. The 
first 12 are for fixed REP. No monotonistically increasing 
behavior was observed for the means and the variances until 
RSR = 0.8. Below RSR of 0.8 the variances are very erratic. 
This is more noticeable in the second 12 graphs where RSR is 
fixed. The distribution means and variances are not linear in 
CRR even when RSR = 1.0. The worst case is when RSR = 0.0. 
Notice that when RSR = 0.0, the largest mean occurs when CRR 
= 0.0, and the smallest occurs when CRR = 0.2. Also, a larger 
mean is achieved when CRR = 0.8 than when CRR = 1.0. In 
addition, the distribution mean is not an increasing function 
of CRR until RSR 2 0.4. The distribution variance has a 
similar behavior to the mean when CRR 2 0.6 and the opposite 
for CRR < 0.6. 

A review of the graphs of the distribution mean as 
a function of CRR when RSR is small suggests an explanation 
which may have some credibility. When CRR = 0, there is only 
one queue affecting the system; i.e., the attrition queue. 
When carcasses are returned (CRR > 0), both the repair and the 
attrition queues affect the system. The drop in the value of 
the mean and the increase in the variance values for RSR < 0.6 
can be attributed to the fact that because the repair queue 


fills at such a slow rate, the delay in identifying non- 
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repairable carcasses whic . then pass to the attrition queue 
allows so much time to pass that demands "eat up" the net 
inventory. In other words, attritions are waiting in the 
repair queue and will not be counted as such until identified 
in the induction process. From an operational point of view 
these results suggest that it is better not to return any 
carcasses when the repair survival rate is very low. 

As the carcass return rate increases for RSR 2 0.4, 
the repair process has a greater effect on the system and, 
because repair turn-around time is shorter than procurement 
lead time, the value of the mean increases. However, RSR of 
at least 0.8 is needed before the values of the variances 
decrease as a function of CRR over its full range of values. 
When CRR = 1.0, again only one queue affects the system; i.e., 
the repair queue. 

b. Distribution Mean and Variance as a function of 
the product of CRR and RSR 

The sensitivity analysis of the effect of the 
product of CRR and RSR on the net inventory’s mean and 
variance was also conducted using the data of Appendix K. 
Because UICP uses (1 - CRR * RSR) to compute the expected 
attrition demand, plots of the distribution mean and variance 
as a function of the product of CRR and RSR were generated for 
fixed values of REP. Because these plots appeared "linear", 


a regression line was fitted to each plot using the AGSS 
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software and, as previously stated, AGSS generated the 
associated ANOVA tables. The plots are provided in Appendix 


N. A summary of the test results are shown in Tables 3 and 4. 





TABLE 3. STATISTICAL SUMMARY RESULTS OF 
THE MEAN VS PRODUCT OF CRR AND RSR. 





a, a DS OE AE SO 
TABLE 4. STATISTICAL SUMMARY RESULTS OF 
THE VARIANCE VS THE PRODUCT OF CRR AND RSR. 


Two tests were performed by the AGSS software, one for the 
regression and the other for the "Lack of Fit". 

A quick look at the fitted lines for the various 
‘values of REP indicate that there does indeed appear to be a 
linear relationship between the distribution mean and the 
product of CRR and RSR. This conclusion is supported, though 


not strongly, by the results of the regression and lack of fit 
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tests. In the regression test the null hypothesis is that the 


distribution mean is not a function of the product; namely, 
the slope of the regression line, B,, is zero. At a = .05 and 
-01, and with the combined degrees of freedom of 1 and 28, the 
critical F-statistics are F,, 1,2. = 4.196 and Foy... = 7.636, 
respectively. Note in Table 3 that in all cases the F- 
statistic (F,), which ranges from 2,608 to 6,690, is much 
greater than these critical values. Therefore, we would 
reject the null hypothesis. 

For the test of lack of fit, the null hypothesis 
becomes H,: E(y) = B, + B,x, a straight line. The alternative 
hypothesis, H,, is that E(y) is some other function of x. At 
a = .05 and .01 with the combined degrees of freedom of 13 and 
15 the critical values are F,, 4315 = 2.45 and Fy,433, = 3.615, 
respectively. Referring to Table 3, we would accept the null 
hypothesis at the .01 significance level because all the F- 
statistics (F,-) are less than the critical value. However, 
for a significance level of .05, we would reject the null 
hypothesis for the cases where REP = 0 and REP = .25 because 
their F-statistics are greater than the critical value. We 
would not reject the null hypothesis in the cases where REP 
= 0.5, 0.75, 1.0 and 1.3 because their F-statistics were less 
than the critical value. The results of these tests indicate 
that the mean is very close to being a linear function of the 
product of CRR and RSR (indeed, the R? values are all very 


close to one), but there is likely at least one other term 
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involved for low REP values. Notice that in all plots the 


CRR * RSR = 1.0 point is the point furthest from the fitted 
line. The reason for this is not obvious. 

The same hypothesis testing as was done for the 
mean was performed for the variance, and the results are 
presented in Table 4. The regression null hypothesis (i.e., 
slope of the regression line is zero) for all values of REP is 
rejected because all the observed F-statistics are greater 
than the critical statistic, F,,,3, = 7.636 for the 0.01 
significance level. The null hypothesis that the variance is 
a linear function of the CRR * RSR product is accepted for all 
values of REP because .the observed F-statistics are all 
greater than the critical F-statistics, Fy 43;,;, = 3.615. 

Finally, the regression lines of the means and 
variances for each value of REP was plotted in the same graph 
to determine if there was any relationship between the mean 
and variance. These graphs are presented in Appendix 0. The 
key relationship indicated from these observations is that the 
variance to mean ration is extremely large when the product of 
CRR and RSR is small; and as the product increases, the ratio 
decreases. Also noted is the fact that the product value 
where the lines intersect (and the variance to mean ratio is 


1.0) increases as REP increases. 
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V. SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


A. SUMMARY 

The current Integrated Repairables Model which is used by 
the Navy’s Inventory Control Points was first conceived in the 
late 1970’s and implemented in 1984. Its development was 
intended as an interim "“patch" to the previous model 
(developed in the late 1960’s) because there tended to be a 
problem of not having enough carcasses available to induct for 
repair when the economic repair quantity and the old repair 
reorder point were used. Its approach was intuitively simple; 
i.e., it provided a common safety stock for both the repair 
and procurement reorder points. This safety stock was derived 
from a RISK formula which was based on a weighted average of 
the unit costs to procure and repair. It also combined the 
procurement lead time and the repair turn-around time with the 
same weights to generate an average lead time for the 
procurement part of the model. A Total Variable Cost equation 
was postulated but the details of all of the terms were never 
worked out so that optimization could justify the formulas for 
the reorder points. In addition, the order quantity formulas 
remained the same as those from the model used prior to the 


integrated model. 





The problem with the old and integrated repairables models 


was that the models as detailed in the Functional Description 
of the CARES (Computation and Research Evaluation System) III 
of UICP were really two models. One model was for procurement 
of new units and the other was for repair of old ones. These 
models had separate reorder points associated with 
"fictitious" inventory positions. The reality is that there 
is only one aggregate inventory made up of new procured and 
old repaired units. 

This thesis takes the first step in modelling that single 
inventory problem. It takes a different approach from the 
classical continuous review inventory model which uses reorder 
points to trigger repair and procurement orders. Instead, it 
collects demands for two “economic order quantities", one for 
the procurement of new items, and the other for the repair of 
the broken ones. It does this by creating two queues; the 
repair eimie which accumulates NRFI carcasses, and the 
attrition queue, which accumulates attritions from either the 
non-receipt of a carcass from a customer or a condemnation of 
a carcass at a maintenance activity. A repair order is 
generated when the repair queue reaches a level corresponding 
to the economic repair quantity, then the contents of that 
queue are reset to zero. The attrition queue behaves 
similarly. Thus, there is no need for a "reorder level" per 


se. 
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The major tasking of this study was to develop a 


simulation program which accurately models the repair process 
and incorporates the queues just mentioned. This simulation 
models the repair process as it really is at an organic depot. 
A carcass is inducted, determined to be repairable or not, and 
then sent either on through repair or sent immediately to the 
attrition queue. Each carcass going through repair takes the 
same specified repair turn-around time to complete the 
process. However, when a batch of carcasses is "inducted", 
each unit must, in fact, wait its turn for induction and that 
wait is not part of its repair turn-around time. 

The purpose of the simulation was to generate the 
appropriate data for use in analyzing the effects of the 
procurement and repair processes on the aggregate inventory 
position probability distribution, the net inventory 
probability distribution, and the safety level. In 
particular, three separate analyses were conducted. In the 
first, inventory position values were simulated for different 
values of the three parameters, the Carcass Return Rate (CRR), 
the Repair Survival Rate (RSR), and the Repair Processing Time 
(REP). The resulting probability distributions were compared 
to the theoretical trapezoidal distribution described in 
Chapter III, which had been developed by Professor A. W. 
McMasters of the Naval Postgraduate School. Next, a way to 
define the aggregate safety level was discovered from studying 


simulation plots of net inventory. Then simulated safety 
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levels were compared to the UICP computed safety levels. 
Finally, net inventory distributions were generated for 
various values of the three previously mentioned parameters as 
a first step in developing an analytical model of the 


aggregate net inventory probability distribution. 


B. CONCLUSIONS. 
1. Inventory Position. 

The results of the comparison of the simulated 
inventory position distributions with the proposed trapezoidal 
distribution were quite revealing. For the limiting cases 
(i.e., when the Carcass Return Rate (CRR) equals zero; CRR 
equals one and the Repair Survival Rate (RSR) is zero; or 
when both CRR and RSR are one.) the trapezoidal model does not 
fit the data. This, however, was expected because under those 
conditions the system parallels the consumable inventory 
system where the inventory position distribution is known to 
be uniformly distributed. [Ref. 3] However, for the "in 
between" cases, (i.e., those where repair and procurement 
processes are both active) the trapezoidal distribution fits 
the data extremely well even though the _ simplifying 
assumptions under which it was derived are not met. Thus, the 
trapezoidal model is quite robust. 

2. Safety Level. 
The results of the comparison of the simulated safety 


levels with those of the UICP model indicate that in all but 
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two cases, the simulated safety level was greater than the 
total safety level (i.e., the sum of the two reorder point 
safety levels of UICP). Thus, the total safety level computed 
in CARES III understates reality. 
3. Net Inventory 

As noted previously, in the Hadley-Whitin consumable 
model which assumes demand to be a Poisson process and lead 
time to be deterministic, one can approximate the demand 
during lead time by using the Normal Distribution with yw equal 
to mean lead time demand and o? = yp. The formula for the net 
inventory distribution is then given in Hadley and Whitin 
(Ref. 3) and it is not Normal. Interestingly though, the 
fitting of the Normal distribution to the various simulated 
net inventory distributions proved to be very good in all 
cases. Although one cannot conclude yet that the net 
inventory distribution can be approximated by the Normal 
distribution for the aggregate repairables model for all 
parameter values affecting the model, this result certainly 
whets the appetite for further investigation. 

From the sensitivity analysis of the three parameters, 
CRR, RSR and REP, it appears that REP has a linear effect on 
the mean and variance of the net inventory distribution; as 
REP increases, the mean decreases, while the variance 
increases most of the time, particularly when CRR and RSR are 


greater than or equal to 0.8. The effects on the mean and 
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variance of the net inventory distribution of the individual 


parameters, CRR and RSR, were not linear. As expected, the 
mean was an increasing function of each parameter except for 
3mall values of CRR. Similarly, the variance was a decreasing 
function of each parameter except for small values of CRR. 
Because the UICP model uses CRR and RSR as a product, a 
sensitivity analysis was conducted to analyze the effect of 
that product on the mean and variance of the net inventory 
distribution. The results showed a statistically stable 
linear form for both, with the mean increasing and the 


variance decreasing as the product value increased. 


C. RECOMMENDATIONS. 

This thesis represents the "tip of the iceberg" in the 
analysis of the aggregate repairables inventory system and the 
development of an analytical model of it. Because the Navy's 
ICP’s intend to become the Repairables "Center of Excellence", 
they must be able to manage such items effectively. The model 
for the theoretical aggregate repairables inventory system can 
aid in that process. 

Now that a simulation model exists, more sensitivity 
analyses of the inventory position and net inventory 
probability distributions must be accomplished. Clearly, the 
effect of low values of CRR on the net inventory distribution 


parameters needs to be better understood. In addition, major 


80 


sensitivity analyses of tnese distributions must include the 


effects caused by the following: 
1. High, medium and low quarterly demand rates; 
2. Actual empirical demand distributions; 
3. Long, medium, and short procurement lead times; 
4. Long, medium, and short repair turn-around times; and, 
5. Stochastic procurement lead times and repair turn-around 


times. 


It is also recommended that the Naval Postgraduate School and 
the Naval Supply Systems Command share in the research effort 


so that it is accomplished in an expeditious manner. 
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APPENDIX A 


MODEL DEFAULTS 


Model Name: THESIS 

Model Description: DEPOT LEVEL REPAIRABLES 
Output File: THSVER4 .OUT 

Run Mode: HIGH_SPEED 

Trace Vars: IP, INV 


Random Number Seed:18645 

Initial Values: 8.2, 1, .8, 1.3, 36, 10, 5, 9, .50 
Ending Condition: STOP_ON_EVENT 

Event: RUN Number to Run: 1 

Trace Events: ALL EVENTS TRACED 


STATE VARIABLES 


wee wr were wrenwe ev 


State Variable #1 


Name: L 

Description: PROCUREMENT LEAD TIME 
Type: REAL 

Size: 1 


State Variable #2 


Name: CRR 

Description: CARCASS RETURN RATE 
Type: REAL 

Size: 1 


State Variable #3 


Name: RSR 

Description: REPAIR SURVIVAL RATE OF A CARCASS. 
Type: REAL 

Size: BS 


State Variable #4 


Name: T2 

Description: REPAIR TURNAROUND TIME 
Type: REAL 

Size: 1 
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State Variable #5 


Name: D 

Description: AVERAGE DEMAND 
Type: REAL 

Size: 1 


State Variable #6 


Name: QP 

Description: PROCUREMENT ORDER QUANTITY 
Type: INTEGER 

Size: 1 


State Variable #7 


Name: QR 

Description: REPAIR ORDER QUANTITY 
Type: INTEGER 

Size: 1 


State Variable #8 


Name: INV 
Description: ACTUAL ON HAND INVENTORY 
Type: INTEGER 
Size 1 
State Variable #9 
Name: IP 
Description: INVENTORY POSITION 
Type: INTEGER 
Size: 1 
State Variable #10 
Name: QUEREP 
Description: NUMBER OF UNITS IN QUEUE FOR REPAIR 
Type: INTEGER 
Size: 1 
State Variable #11 
Name: QUEORD 
Description: NUMBER OF UNITS ON ORDER 
Type: INTEGER 
Size: 1 
State Variable #12 
Name: RQUE 
Description: NUMBER OF UNITS IN QUEUE FOR A REPAIR ORDER 
Type: INTEGER 
Size: 1 
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State Variable #13 


Name: TIME 

Description: LENGTH OF TIME FOR ONE RUN OF THE SIMULATION. 
Type: REAL 

Size: 1 


State Variable #14 


Name: RUN 
Description: COUNTER FOR THE NUMBER OF RUNS 
Type: REAL 
Size: 1 

State Variable #15 
Name: ATTR 
Description: CUMULATIVE TOTAL OF CONDEMN/LOST CARCASSES 
Type: INTEGER 
Size: 1,1 

State Variable #16 
Name: Xx 
Description: HOLDER OF A RANDOM VARIABLE FOR A STATE. 
Type: REAL 
Size: 3 

State Variable #17 
Name: REP 
Description: EVALUATION PROCESSING TIME FOR REPAIR. 
Type: REAL 
Size: 1,1 

State Variable #18 
Name: : COUNT 
Description: NUMBER OF REPAIRABLE CARCASSES IN A BATCH 
Type : INTEGER 
Size: 5 ee 

VERTICES 

Vertex #1 

Name: RUN 


Description: START OF THE SIMULATION 
State Changes: RUN=1 
Parameter(s): L,CRR,RSR,T2,TIME,QP,QR,D,REP 
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Vertex #2 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #3 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #4 
Name: 
Description: 
State Changes: 
Parameter (s): 


Vertex #5 
Name: 
Description: 
State Changes: 


Parameter (s) : 


Vertex #6 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #7 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #8 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #9 
Name: 
Description: 
State Changes: 
Parameter (s) : 





DMAND 
DEMAND OF A DEPOT LEVEL REPAIRABLE 
INV=INV-1,X[1] =RND, IP=IP-1 


ORDER 
ORDERING OF QP WNITS OF NEW MATERIAL 
IP=IP+QP, QUEORD=QUEORD+QOQP, ATTR=0 


RCVD2 
RECEIPT OF AN ORDER OF QP UNITS 
INV=INV+QP, QUEORD=QUEORD-QP 


CARTN 

CARCASS RETURNING TO THE STOCK POINT 

RQUE=RQUE+1, COUNT=1, IP=IP+RQUE* (RQUE==QR) , 
QUEREP=QUEREP+ (RQUE==QR) *RQUE, X [2] =0 


RPORD 
DECISION TO REPAIR OR CONDEMN A CARCASS 
X [0] =RND, RQUE=RQUE-1 


DEPOT 
DEPOT REPAIRING A CARCASS 


ATTRI 
LOSS OR CONDEMNATION OF A CARCASS 
ATTR«ATTR+1 


STATS 
ACCUMULATION OF SIMULATION RUNS 
RUN=RUN+1 
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Vertex #10 
Name: INIT 
Description: INITIALIZE iP AND INV FOR EACH RUN 
State Changes: INV=72,IP=INV 
Parameter (s) : 


Vertex #11 
Name: ATTR1 
Description: CAPTURING OF THE IP JUST BEFORE AN ATTRITION 
State Changes: QUEREP=QUEREP-1,X[2] =1 
Parameter (s): 


Vertex #12 
Name : RCVD1 
Description: RECORDING OF THE INV VALUE JUST BEFORE A 
RECEIPT 
State Changes: 
Parameter (s): 
Vertex #13 
Name: SHIP1 


Description: RECEIPT OF A REPAIRED CARCASS 
State Changes: INV=INV+1. 
Parameter (s): 


Vertex #14 
Name: DMND1 
Description: CAPTURING OF THE INV AND IP JUST BEFORE A 
DEMAND 


State Changes: 
Parameter (s) : 


Vertex #15 
Name: RPAIR 
Description: INV VALUE JUST BEFORE RECEIPT OF A REPAIRED 
UNIT 
State Changes: QUEREP=QUEREP-1 
Parameter (s) : 


Vertex #16 
Name : ATTR2 
Description: REDUCTION OF IP DUE TO A REPAIR ATTRITION 
State Changes: ATTReATTR+1,IP=IP-1 
Parameter (8s) : 


Vertex #17 
Name: COUNT 
Description: NUMBER OF SUCCESSFUL REPAIRS IN A BATCH 
State Changes: COUNT=COUNT+ (X [0] <=RSR) 
Parameter (s) : 
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Vertex #18 
Name : INPOS 
Description: INCREASE OF IP BY INDUCTING A BATCH IN REPAIR 
State Changes: IP=IP 
Parameter (s) : 


EDGES 
Edge #1 

Description: A CARCASS IS SUCCESSFULLY RETURNED 
Type: Scheduling 
Origin: DMAND 
Destination: CARTN 
Condition: X [1] <=CRR 
Delay: 0 
Priority: 3 
Attributes: 

Edge #2 
Description: BEGIN THE SIMULATION 
Type: Scheduling 
Origin: RUN 
Destination: INIT 
Condition: l==1 
Delay: 0 
Priority: 5 
Attributes: 

Edge #3 
Description: START THE NEXT SIMULATION 
Type: Scheduling 
Origin: STATS 
Destination: INIT 
Condition: l==1 
Delay: fe) 
Priority: 5 
Attributes: 

Edge #4 
Description: END THE RUN OF THE SIMULATION 
Type: Scheduling 
Origin: INIT 
Destination: STATS 
Condition: 1l==1 
Delay: TIME 
Priority: 5 
Attributes: 
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Edge #5 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #6 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #7 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #8 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #9 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


CANCEL ALL FUTURE DEMANDS 


Cancelling 
STATS 


A CARCASS IS NOT RETURNED TO THE STOCK POINT 


Scheduling 
DMAND 


PROCURE QP UNITS 


Scheduling 
ATTRI 
ORDER 
ATTR>=QP 

0 

3 


SCHEDULE THE RECEIPT OF QP UNITS 


Scheduling 
ORDER 
RCVD1 

le=1 

L 

5 


ADD QP TO INVENTORY (INV) 


Scheduling 
RCVD1 
RCVD2 

le=«1 

0 

5 
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Edge #10 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #11 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #12 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #13 
Description:. 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #14 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Cancelling 
STATS 
RCVD1 

issl 

0 

5 


* 





BEGIN THE DEMAND CYCLE OF THE SIMULATION 
Scheduling 

INIT 

DMND1 

lesl 

0 

5 


SCHEDULE ANOTHER DEMAND 
Scheduling 


l==1 
- (1/D) *LN{RND} 
5 


RECORD INV AND IP JUST PRIOR TO A DEMAND 
Scheduling 
DMND1 


SHIP THE RFI CARCASS TO THE STOCK POINT 
Scheduling : 
DEPOT 
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Edge #15 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #16 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #17 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #18 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #19 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


RECEIVE A JUST REPAIRED UNIT 
Scheduling 

RPAIR 

SHIP1 

i=sl1 

0 

5 


SCHEDULE THE REPAIR ATTRITION 
Scheduling 
RPORD 
ATTR1 
X [0] >RSR 
REP* COUNT 
3 


SCHEDULE THE REPAIR OF A NRFI UNIT 
Scheduling 
RPORD 
DEPOT 
X [0] <=RSR 
REP* COUNT 
3 


SEND A BATCH OF QP NRFI CARCASSES TO A DEPOT 
Scheduling 

CARTN 

RPORD 

RQUE>=QR 

0 

5 


INCREASE ATTR AND REDUCE IP BY ONE 
Scheduling 
ATTR1 
ATTR2 
lel 
0 
7 
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Edge #20 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #21 
Description: 
BATCH 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #22 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #23 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


INCREASE THE COUNT OF RFI UNITS IN THE BATCH 
Scheduling 
RPORD 


SCHEDULE REVIEW OF THE NEXT CARCASS IN THE 
Scheduling 

COUNT 

RPORD 

RQUE>0 

0 

5 


PROCURE QP UNITS 
Scheduling 
ATTR2 
ORDER 
ATTR>=QP 
0 
2 


RECORD THE IP AS A RESULT OF THIS BATCH 
Scheduling 
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APPENDIX B 


The SIGMA Model, THESIS, is a discrete event simulation. It 
models depot level repairables. 


I. STATE VARIABLE DEFINITIONS. 


For this simulation, the following state variables are 
defined: 

L: procurement lead time (real valued) 

CRR: carcass return rate (real valued) 

RSR: repair survival rate of a carcass. (real valued) 

T2: repair turnaround time (real valued) 

D: average demand (real valued) 

QP: procurement order quantity (integer valued) 

QR: repair order quantity (integer valued) 

INV: actual on hand inventory (integer valued) 

IP: inventory position (integer valued) 

QUEREP: number of units in queue for repair (integer 
valued) 

QUEORD: number of units on order (integer valued) 

RQUE: number of units in queue for a repair order 
(integer valued) 

TIME: length of time for one run of the simulation. (real 

valued) 
RUN: counter for the number of runs’ (real valued) 
ATTR: cumulative total of condemn/lost carcasses (integer 
valued) 

X[(3]: holder of a random variable for a state. (real 
valued) 

REP: evaluation processing time for repair. (real valued) 

COUNT: number of repairable carcasses in a batch 
(integer valued) 


II. EVENT DEFINITIONS. 


Simulation state changes are represented by event vertices 
(nodes or balls) in a SIGMA graph. Event vertex parameters, 
if any, are given in parentheses. 

Logical and dynamic relationships between pairs of events 
are represented in a SIGMA graph by edges (arrows) between 
event vertices. Unless otherwise stated, vertex execution 
priorities, to break time ties, are equal to 5. 
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1. The RUN(L,CRR,RSR,T2,TIME,QP,QR,D,REP) event models the 
start of the simulation. 
Initial values for, L, CRR, RSR, T2, TIME, QP, QR, D, REP, 
are needed for each run. 
This event causes the following state change(s) : 
RUN=1 
After every occurrence of the RUN event: 
Unconditionally, begin the simulation; 
that is, schedule the INIT() event to occur without 
delay. 


2. The DMAND() event models the demand of a depot level 
repairable. 

This event causes the following state change(s) : 
INV=INV-1 
X (1) =RND 
IPxIP-1 

After every occurrence of the DMAND event: 

If X[(1]<=CRR, then a carcass is successfully returned; 
that is, schedule the CARTN() event to occur without 
delay. 

(Time ties are broken by an execution priority of 3.) 

If CRR<X[1], then a carcass is not returned to the stock 

point; 
that is, schedule the ATTRI() event to occur without 
delay. 
(Time ties are broken by an execution priority of 3.) 

Unconditionally, schedule another demand; 
that is, schedule the DMND1() event to occur in 

-(1/D)*LN{RND} time units. 


3. The ORDER() event models the ordering of QP units of new 
material. 

This event causes the following state change(s): 
IP=IP+QP 
QUEORD=QUEORD+QP 
ATTR=0 

After every occurrence of the ORDER event: 

Unconditionally, schedule the receipt of QP units ; 
that is, schedule the RCVD1() event to occur in L 
time units. 


4. The RCVD2() event models the receipt of an order of QP 
units. 
This event causes the following state change(s): 
INV=INV+QP 
QUEORD=QUEORD-QP 
No additional events are scheduled here. 
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5. The CARTN() event models the carcass returning to the stock 
point. 

This event causes the following state change(s): 
RQUE=#RQUE+1 
COUNT=1 
IP=2IP+RQUE* (RQUE==QR) 

QUEREP=QUEREP+ (RQUE==QR) *RQUE 
X[2] =0 
After every occurrence of the CARTN event: 
If RQUE>=QR, then send a batch of QP NRFI carcasses to a 
depot ; 
that is, schedule the RPORD() event to occur without 
delay. 

If RQUE==<QR, then record the IP as a result of this batch; 
that is, schedule the INPOS() event to occur without 
delay. 

(Time ties are broken by an execution priority of 8.) 


6. The RPORD() event models the decision to repair or condemn 
a carcass. 
This event causes the following state change(s): 
X [0] =RND 
RQUE=RQUE-1 
After every occurrence of the RPORD event: 
If X[0]>RSR, then schedule the repair attrition ; 
that is, schedule the ATTR1() event to occur in 
REP* (COUNT-1) 
time units. 
(Time ties are broken by an execution priority of 3.) 
If X[0)<=#RSR, then schedule the repair of a NRFI unit; 
that is, schedule the DEPOT() event to occur in 
REP* (COUNT-1) 
time units. 
(Time ties are broken by an execution priority of 3.) 
Unconditionally, increase the count of RFI units in the 
batch; 
that is, schedule the COUNT() event to occur without 
delay. 


7. The DEPOT() event models the depot repairing a carcass. 
After every occurrence of the DEPOT event: 
Unconditionally, ship the RFI carcass to the stock point; 

that is, schedule the RPAIR() event to occur in T2 
time units. 
(Time ties are broken by an execution priority of 3.) 





8. The ATTRI() event models the loss or condemnation of a 
carcass. 
This event causes the following state change(s): 
ATTR=ATTR+1 
After every occurrence of the ATTRI event: 
If ATTR>=#QP, then procure QP units; 
that is, schedule the ORDER() event to occur without 
delay. 
(Time ties are broken by an execution priority of 3.) 


9. The STATS() event models the accumulation of simulation 
runs. 
This event causes the following state change(s) : 
RUN=RUN+1 
After every occurrence of the STATS event: 
Unconditionally, start the next simulation; 
that is, schedule the INIT() event to occur without 
delay. 
Unconditionally, cancel all future demands; 
that is, immediately cancel all scheduled occurrences 
of the DMAND event. 
Unconditionally, immediately cancel all scheduled 
occurrences 
of the RCVD1 event. 





10. The INIT() event models the initialize IP and INV for each 
run. 
This event causes the following state change(s): 
INV=72 
IP=INV 
After every occurrence of the INIT event: 
Unconditionally, end the run of the simulation.; 
that is, schedule the STATS() event to occur in TIME 
time units. 
Unconditionally, begin the demand cycle of the simulation; 
that is, schedule the DMND1() event to occur without 
delay. 


11. The ATTR1() event models the capturing of the IP just 
before an attrition. 
This event causes the following state change(s) : 
QUEREP=QUEREP-1 
X[(2)=1 
After every occurrence of the ATTR1 event: 
Unconditionally, increase ATTR and reduce IP by one; 
that is, schedule the ATTR2() event to occur without 
delay. 
(Time ties are broken by an execution priority of 7.) 
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12. The RCVD1() event models the recording of the INV value 
just before a receipt. 
After every occurrence of the RCVD1 event: 
Unconditionally, add QP to inventory (INV); 
that is, schedule the RCVD2() event to occur without 
delay. 


13. The SHIP1() event models the receipt of a repaired 
carcass. 
This event causes the following state change(s): 
INV=INV+1 
No additional events are scheduled here. 


14. The DMND1() event models the capturing of the INV and IP 
just before a demand. 
After every occurrence of the DMND1 event: 
Unconditionally, record INV and IP just prior to a demand; 
that is, schedule the DMAND() event to occur without 
delay. 


15. The RPAIR() event models the INV value just before receipt 
of a repaired unit. 
This event causes the following state change(s) : 
QUEREP=QUEREP-1 
After every occurrence of the RPAIR event: 
Unconditionally, receive a just repaired unit; 
that is, schedule the SHIP1() event to occur without 
delay. 


16. The ATTR2() event models the reduction of IP due to a 
repair attrition. 
This event causes the following state change(s): 
ATTR=ATTR+1 
IP=IP-1 
After every occurrence of the ATTR2 event: 
If ATTR>=QP, then procure QP units; 
that is, schedule the ORDER() event to occur without 
delay. 
(Time ties are broken by an execution priority of 2.) 


17. The COUNT() event models the number of successful repairs 
in a batch. 
This event causes the following state change(s) : 
COUNT=COUNT+ (X [0] <=RSR) 
After every occurrence of the COUNT event: 
If RQUE>0, then schedule review of the next carcass in the 
batch; 
that is, schedule the RPORD() event to occur without 
delay. 
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18. The INPOS() event models the increase of IP by inducting 
a batch in repair. 
This event causes the following state change(s): 
TPp=IP 
No additional events are scheduled here. 


97 











APPENDIX C 


MODEL DEFAULTS 


Model Name: THESIS 

Model Description: DEPOT LEVEL REPAIRABLES 
Output File: THSVER4 .OUT 

Run Mode: HIGH_SPEED 

Trace Vars: IP, INV 

Random Number Seed:18645 

Initial Values: 8.2, .6, .6, 1.3,36, 


Ending Condition: STOP_ON_EVENT 
Event: STATS Number to Run: 1 


Trace Events: ALL EVENTS TRACED 

Time Event Count Ip INV 
0.000 RUN 1 0 0 
0.000 INIT 2 72 72 
0.000 DMND1 1 72 72 
0.000 DMAND 1 71 71 
0.000 CARTN 1 71 71 
0.072 DMND1 2 71 71 
0.072 DMAND 2 70 70 
0.072 CARTN 2 70 70 
0.142 DMND1 3 70 70 
0.142 DMAND 3 69 69 
0.142 CARTN 3 69 69 
0.221 DMNDi 4 69 69 
0.221  DMAND 4 68 68 
0.221 CARTN 4 68 66 
0.314 DMND1 5 68 68 
0.314 DMAND 5 67 67 
0.314 CARTN 5 72 67 
0.314 RPORD 1 72 67 
0.314 DEPOT 2 72 67 
0.314 COUNT 1 72 67 
0.314 RPORD 2 72 67 
0.314 COUNT 2 72 67 
0.314 RPORD 3 72 67 
0.314 COUNT 3 72 67 
0.314 RPORD 4 72 67 
0.314 COUNT 4 72 67 
0.314 RPORD 5 72 67 
0.314 COUNT 5 72 67 
0.314 INPOS 1 72 67 
0.431 DMNDi 6 72 67 
0.431 DMAND 6 71 66 
0.431 CARIN 6 71 66 
0.491 DMND1 7 71 66 
0.491 DMAND 7 70 65 
0.491 CARTN 7 70 65 
0.518 DMND1 8 70 65 
0.518 DMAND 8 69 64 
0.518 CARTN 8 69 64 
98 
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1.569 
1.569 
2.569 
2.577 
1.577 
1.614 
1.614 
1.733 
1.733 
1.733 
1.819 
1.619 
1.619 
1.819 
1.819 
1.819 
1.864 
1.864 
1.865 
1.665 
1.865 
1.943 
1.943 
1.943 
2.114 
2.114 
2.182 
2.182 
2.182 
2.364 
2.364 
2.369 
2.369 
2.369 
2.377 
2.377 
2.520 
2.520 
2.520 
2.627 
2.627 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.668 
2.675 
2.675 
2.675 
2.798 
2.798 
2.798 
2.869 
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APPENDIX D 


MODEL DEFAULTS 


Model Name: THESIS 

Model Description: DEPOT LEVEL REPAIRABLES 
Output File: THSVER7 .OUT 

Run Mode: HIGH_SPEED 

Trace Vars: CRR,RSR, REP, IP 

Random Number Seed:18645 

Initial Values: 8.2, 1.3, 850, 10, 5S, 10 


Ending Condition: STOP_ON_EVENT 
Event: STATS Number to Run: 1 
Trace Events: ALL EVENTS TRACED 


STATE VARIABLES 


were wr ewww weenece 


State Variable #1 


Name: L 

Description: PROCUREMENT LEAD TIME 
Type: REAL 

Size: 1 


State Variable #2 


Name: CRR 

Description: CARCASS RETURN RATE 
Type: REAL 

Size: 1 


State Variable #3 


Name: RSR 

Description: REPAIR SURVIVAL RATE OF A CARCASS. 
Type: REAL 

Size: 1 


State Variable #4 


Name: T2 

Description: REPAIR TURNAROUND TIME 
Type: REAL 

Size: 1 
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State Variable #5 


Name : D 

Description: AVERAGE QUARTERLY DEMAND 
Type: REAL 

Size: 1 


State Variable #6 


Name: QP 

Description: PROCUREMENT ORDER QUANTITY 
Type: INTEGER 

Size: ni 


State Variable #7 


Name : QR 
Description: REPAIR ORDER QUANTITY 
Type: INTEGER 
Size: 1 
State Variable #8 
Name: INV 
Description: ACTUAL ON HAND INVENTORY 
Type: INTEGER 
Size: 1 
State Variable #9 
Name: IP 
Description: INVENTORY POSITION 
Type: INTEGER 
Size: 1 


State Variable #10 


Name: QUEREP 

Description: NUMBER OF UNITS IN QUEUE FOR REPAIR 
Type: INTEGER 

Size: 1 


State Variable #11 


Name: QUEORD 

Description: NUMEER OF UNITS ON ORDER 
Type: INTEGER 

Size: 1 


State Variable #12 


Name: RQUE 

Description: NUMBER OF CARCASSES AWAITING REPAIR ORDER 
Type: INTEGER 

Size: 1 
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State Variable #13 


Name : TIME 

Description: LENGTH OF TIME FOR ONE RUN OF THE SIMULATION. 
Type: REAL 

Size: 1 


State Variable #14 


Name: SAFTY 

Description: COUNTER FOR THE SAFETY LEVELS 
Type: REAL 

Size: 7 


State Variable #15 


Name: ATTR 
Description: CUMULATIVE TOTAL OF CONDEMN/LOST CARCASSES 
Type: INTEGER 
Size: 1,1 
State Variable #16 
Name: x 
Description: HOLDER OF A RANDOM VARIABLE FOR A STATE. 
Type: REAL 
Size: 3 


State Variable #17 


Name: REP 
Description: EVALUATION PROCESSING TIME FOR REPAIR. 
Type: REAL 
Size: 1,1 

State Variable #18 
Name: COUNT 
Description: COUNTER FOR THE REPAIR PROCESS 
Type: INTEGER 
Size: 1,1 

State Variable #19 
Name: COUNT1 
Description: COUNTER FOR THE NUMBER OF INV CHANGES 
Type: INTEGER 
Size: 1 

VERTICES 

Vertex #1 

Name: RUN 


Description: START OF THE SIMULATION 
State Changes: 
Parameter (s): L,T2, TIME, QP,QR 
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Vertex #2 
Name: 
Description: 
State Changes: 
Parameter (3) : 


Vertex #3 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #4 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #5 
Name : 
Description: 
State Changes: 


Parameter (s) : 


Vertex #6 
Name: 
Description: 
State Changes: 
Parameter (s): 


Vertex #7 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #8 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #9 
Name: 
Description: 
State Changes: 
Parameter (s) : 


DMAND 
DEMAND OF A DEPOT LEVEL REPAIRABLE 
INV=INV-1,X[1] =RND, IP=IP-1 


ORDER 
PROCURING OF QP UNITS 
IP=IP+QP, QUEORD=QUEORD+QP, ATTR=0 


RCVD2 
RECEIPT OF AN ORDER OF QP UNITS 
INV=INV+QP, QUEORD=QUEORD-QP 


CARTN 

CARCASS RETURNING TO THE STOCK POINT 

RQUE=RQUE+1, COUNT=0, IP=IP+RQUE* (RQUE==QR) , 
QUEREP=QUEREP+ (RQUE==<QR) *RQUE, X[2]=0 


RPORD 
DECISION TO REPAIR OR CONDEMN A CARCASS 
X [0] =RND, RQUE=RQUE-1 


DEPOT 
DEPOT REPAIRING A CARCASS. 


ATTRI 
LOSS OR CONDEMNATION OF A CARCASS 
ATTR=ATTR+1 


STATS 
COMPUTATION OF AVERAGE SAFETY LEVEL 
SAFTY=SAFTY/COUNT1 
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Vertex #10 
Name : INIT 
Description: INITIALIZATION OF IP, INV, D, CRR, RSR, AND REP 


s t a t e Cc h a n g e 
INV=72, IP=INV, D=DISK{THS.DAT;0},CRR=DISK{THS.DAT;0}, 

RSR=DISK{THS.DAT;0}, REP=DISK{THS.DAT;0} 
Parameter (s) : 


s 


Vertex #11 
Name: ATTR1 
Description: CAPTURING OF IP JUST BEFORE A REPAIR ATTRITION 
State Changes: QUEREP=QUEREP-1,X[2]=1 
Parameter (s) : 


Vertex #12 
Name: RCVD1 
Description: ACCUMULATION OF INV JUST BEFORE A RECEIPT 
Ss t a t e Cc h a n g e s 


SAFTY=SAFTY+ (CLK>=20) * INV, COUNT1=COUNT1+ (CLK>=20.0) 
Parameter (s) : 


Vertex #13 
Name: SHIP1 
Description: RECEIPT OF A REPAIRED CARCASS 
State Changes: INV=INV+1 
Parameter (s) : 


Vertex #14 
Name: DMND1 
Description: CAPTURING OF INV AND IP JUST PRIOR TO A DEMAND 
State Changes: 
Parameter (s) : 


Vertex #15 


Name: RPAIR 
Description: INV VALUE JUST BEFORE RECEIPT OF A REPAIRED 
UNIT 


State Changes: QUEREP=QUEREP-1,SAFTY=SAFTY+ (CLK>=20) * INV, 
COUNT1=COUNT1+ (CLK>=20) 
Parameter (s) : 


Vertex #16 
Name: ATTR2 
Description: ATTRITION DUE TO REPAIR 
State Changes: ATTR=ATTR+1,IP=IP-1 
Parameter (s) : 
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Vertex #17 


Name: 
Description: 


COUNT 


NUMBER OF SUCCESSFUL REPAIRS IN A BATCH 


State Changes: COUNT=COUNT+ (X [0] <=RSR) 


Parameter (s) : 


Vertex #18 


Name : 
Description: 


State Changes: 
Parameter (s): 


INPOS 


INDUCTION OF A BATCH OF CARCASSES INTO REPAIR 


IP=IP 


Vertex #19 


Name: 
Description: 


Edge #1 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #2 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #3 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


RINIT 


REINITIALIZATION OF TIME, L, T2, QP, AND QR 
State Changes: TIME=850,L=8.2,T2=1.3,QP=10, QR=5 
Parameter (s) : 


A CARCASS IS SUCCESSFULLY RETURNED 


Scheduling 
DMAND 
CARTN 

X [1] <=CRR 
0 

3 


BEGIN THE SIMULATION 
Scheduling 


END THE RUN OF THE SIMULATION 
Scheduling 

INIT 

STATS 

TIME 

4 
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Edge #4 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #5 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #6 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #7 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #8 
Description: 
Type: 
Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 





A CARCASS IS NOT RETURNED TO THE STOCK POINT 
Scheduling 
DMAND 
ATTRI 
CRR<X [1] 
0 
3 


ORDER AN ECONOMIC ORDER QUANTITY. 
Scheduling 
ATTRI 
ORDER 
ATTR>=QP 
0 
3 


MARK THE INVENTORY QTY JUST BEFORE A RECEIPT. 
Scheduling 

ORDER 

RCVD1 

l== 

L 

5 


ADD AN ORDER QTY TO NET INVENTORY 
Scheduling 


SCHEDULE THE NEXT DEMAND 
Scheduling 
DMAND 
DMND1 
1l== 
- (1/D) *LN{RND} 
5 
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Edge #9 


Description: RECORD INV AND IP JUST PRIOR TO DEMAND 
Type: Scheduling 

Origin: DMND1 

Destination: DMAND 

Condition: == 

Delay: 0 

Priority: 5 

Attributes: 

Edge #10 
Description: SHIP THE RFI CARCASS TO THE STOCK POINT. 
Type: Scheduling 
Origin: DEPOT 
Destination: RPAIR 
Condition: == 
Delay: T2 
Priority: 3 
Attributes: 

Edge #11 
Description: INCREASE ATTR AND DECREASE QUEREP BY ONE 
Type: Scheduling 
Origin: RPORD 
Destination: ATTR1 
Condition: X[0] >RSR 
Delay: REP*COUNT 
Priority: 3 
Attributes: 

Edge #12 
Description: 

Type: Scheduling 
Origin: RPORD 
Destination: DEPOT 
Condition: X [0] <=RSR 
Delay: REP* COUNT 
Priority: 3 
Attributes: 

Edge #13 
Description: SUPPLY CENTER SENDS ALL ITS CARCASSES TO 
DEPOT. 

Type: Scheduling 
Origin: CARTN 
Destination: RPORD 
Condition: RQUE>=QR 
Delay: 0 
Priority: 5 
Attributes: 
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Edge #14 


Description: INCREASE ATTR AND REDUCE IP BY ONE 
Type: Scheduling 
Origin: ATTR1 
Destination: ATTR2 
Condition: l==1 
Delay: 0 
Priority: 7 
Attributes: 

Edge #15 
Description: 
Type: Scheduling 
Origin: RPORD 
Destination: COUNT 
Condition: le= 
Delay: 0 
Priority: 5 
Attributes: 

Edge #16 
Description: 
Type: Scheduling 
Origin: COUNT 
Destination: RPORD 
Condition: RQUE>0 
Delay: 0 
Priority: 5 
Attributes: 

Edge #17 
Description: 
Type: Scheduling 
Origin: COUNT 
Destination: INPOS 
Condition: RQUE==0&X [2] == 
Delay: 0 
Priority: 8 
Attributes: 

Edge #18 
Description: 
Type: Scheduling 
Origin: ATTR2 
Destination: ORDER 
Condition: ATTR>=QP 
Delay: 0 
Priority: 2 
Attributes: 
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Edge #19 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #20 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #21 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #22 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 





Scheduling 
RPAIR 
SHIP1 

l==1 

0 

5 


RESET THE SIMULATION MODEL TO ZERO 
Scheduling 

STATS 

RINIT 

1==1&SET{18645} 

0 

8 


BEGIN THE NEXT RUN OF THE SIMULATION 
Scheduling 

RINIT 

INIT 

1==1 

0 

5 


BEGIN THE DEMAND CYCLE 
Scheduling 

INIT 

DMND1 

l== 

0 

5 





APPENDIX E 


The SIGMA Model, THESIS, is a discrete event simulation. It 
models depot level repairables. 


I. STATE VARIABLE DEFINITIONS. 


For this simulation, the following state variables are 
defined: 
L: procurement lead time (real valued) 
CRR: carcass return rate (real valued) 
RSR: repair survival rate of a carcass. (real valued) 
T2: repair turnaround time (real valued) 
D: average quarterly demand (real valued) 


QP: procurement order quantity (integer valued) 

QR: repair order quantity (integer valued) 

INV: actual on hand inventory (integer valued) 

IP: inventory position (integer valued) 

QUEREP: number of units in queue for repair (integer 
valued) 

QUEORD: number of units on order (integer valued) 


RQUE: number of carcasses awaiting repair order 
(integer valued) 

TIME: length of time for one run of the simulation. 
(real valued) 

SAFTY: counter for the safety levels (real valued) 

ATTR: cumulative total of condemn/lost carcasses 
(integer valued) 


X[3]: holder of a random variable for a state. (real 
valued) 

REP: evaluation processing time for repair. (real valued) 

COUNT: counter for the repair process (integer valued) 

COUNT1: counter for the number of inv changes (integer 
valued) 


II. EVENT DEFINITIONS. 


Simulation state changes are represented by event vertices 
(nodes or balls) in a SIGMA graph. Event vertex parameters, 
if any, are given in parentheses. 

Logical and dynamic relationships between pairs of events 
are represented in a SIGMA graph by edges (arrows) between 
event vertices. Unless otherwise stated, vertex execution 
priorities, to break time ties, are equal to 5. 
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1. The RUN(L,T2,TIME,QP,QR) event models the start of the 
simulation. 
Initial values for, L,T2,TIME,QP,QR, are needed for each 
run. 
After every occurrence of the RUN event: 
Unconditionally, begin the simulation; 
that is, schedule the INIT() event to occur without 
delay. 


2. The DMAND() event models the demand of a depot level 
repairable. 

This event causes the following state change(s): 
INV=INV-1 
X(1] =RND 
IP=IP-1 

After every occurrence of the DMAND event: 

If X{1]<=CRR, then a carcass is successfully returned; 
that is, schedule the CARTN() event to occur without 
delay. 

(Time ties are broken by an execution priority of 3.) 

If CRR<X[(1], then a carcass is not returned to the stock 

point; 
that is, schedule the ATTRI() event to occur without 
delay. 
(Time ties are broken by an execution priority of 3.) 

Unconditionally, schedule the next demand; 
that is, schedule the DMND1() event to occur in 
-(1/D) *LN{RND} time units. 


3. The ORDER() event models the procuring of QP units. 
This event causes the following state change(s) : 
IP=IP+QP 
QUEORD=QUEORD+QP 
ATTR=0 
After every occurrence of the ORDER event: 
Uncenditionally, mark the inventory qty just before a 
receip® ; 
tuat is, schedule the RCVD1() event to occur in L 
time units. 


4. The RCVD2() event models the receipt of an order of QP 
units. 
This event causes the following state change(s): 
INV=INV+QP 
QUEORD=QUEORD-QP 
No additional events are scheduled here. 
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5. The CARTN() event models the carcass returning to the stock 
point. 
This event causes the following state change(s) : 
RQUE=RQUE+1 
COUNT=0 
IP=IP+RQUE* (RQUE==QR) 
QUEREP=QUEREP+ (RQUE==QR) *RQUE 
X[2] =0 
After every occurrence of the CARTN event: 
If RQUE>=QR, then supply center sends all its carcasses to 
depot; that is, schedule the RPORD() event to occur without 
delay. 


6. The RPORD() event models the decision to repair or condemn 
a carcass. 
This event causes the following state change(s): 
X [0] =RND 
RQUE=RQUE-1 
After every occurrence of the RPORD event: 
If X[0] >RSR, then increase ATTR and decrease QUEREP by one; 
that is, schedule the ATTR1() event to occur in 
REP* COUNT 
time units. 
(Time ties are broken by an execution priority of 3.) 
If X[0]<=RSR, then schedule the DEPOT() event to occur in 
REP*COUNT time units. 
(Time ties are broken by an execution priority of 3.) 
Unconditionally, schedule the COUNT() event to occur 
without delay. 


7. The DEPOT() event models the depot repairing a carcass. 
After every occurrence of the DEPOT event: 
Unconditionally, ship the RFI carcass to the stock point; 

that is, schedule the RPAIR() event to occur in T2 
time units. 
(Time ties are broken by an execution priority of 3.) 


8. The ATTRI() event models the loss or condemnation of a 
carcass. 
This event causes the following state change(s): 
ATTR=ATTR+1 ; 
After every occurrence of the ATTRI event: 
If ATTR>=QP, then order QP units; 
that is, schedule the ORDER() event to occur without 
delay. 
(Time ties are broken by an execution priority of 3.) 
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9. The STATS() event models the computation of average safety 
level. 
This event causes the following state change(s) : 
SAFTY=SAFTY/COUNT1 
After every occurrence of the STATS event: 
If l==1 and SET{18645}, then reset the simulation model to 
zero; 
that is, schedule the RINIT() event to occur without 
delay. 
(Time ties are broken by an execution priority of 8.) 


10. The INIT() event models the initialization of IP, INV, D, 
CRR, RSR and REP. 
This event causes the following state change(s): 
INV=72 
IP=INV 
D=DISK{THS.DAT; 0} 
CRR=DISK{THS.DAT;0 
RSR=DISK{ THS.DAT; 0 
REP=DISK{THS.DAT;0 
After every occurrence of the INIT event: 
Unconditionally, end the run of the simulation; 
that is, schedule the STATS() event to occur in TIME 
time units. 
(Time ties are broken by an execution priority of 4.) 
Unconditionally, begin the demand cycle; 
that is, schedule the DMND1() event to occur without 
delay. 








11. The ATTR1() event models the capturing of IP just before 
a repair attrition. 
This event causes the following state change(s) : 
QUEREP=QUEREP-1 
X[(2J=1 
After every occurrence of the ATTR1 event: 
Unconditionally, increase ATTR and reduce IP by one; 
that is, schedule the ATTR2() event to occur without 
delay. 
(Time ties are broken by an execution priority of 7.) 


12. The RCVD1() event models the accumulation of INV just 
before a receipt. 
This event causes the following state change(s) : 
SAFTY=SAFTY+ (CLK>=20) * INV 
COUNT1=COUNT1+ (CLK>=20.0) 
After every occurrence of the RCVD1 event: 
Unconditionally, add an order gty to net inventory; 
that is, schedule the RCVD2() event to occur without 
delay. 
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13. The SHIP1() event models the receipt of a repaired 
carcass. 
This event causes the following state change(s) : 
INV=INV+1 
No additional events are scheduled here. 


14. The DMND1() event models the capturing of INV and IP just 
prior to a demand. 
After every occurrence of the DMND1 event: 
Unconditionally, record INV and IP just prior to demand; 
that is, schedule the DMAND() event to occur without 
delay. 


15. The RPAIR() event models the INV value just before receipt 
of a repaired unit. 
This event causes the following state change(s): 
QUEREP=QUEREP-1 
SAFTY=SAFTY+ (CLK>=20) * INV 
COUNT 1=COUNT1+ (CLK>=20) 
After every occurrence of the RPAIR event: 
Unconditionally, schedule the SHIP1() event to occur 
without delay. 


16. The ATTR2() event models the attrition due to repair. 
This event causes the following state change(s) : 
ATTR=ATTR+1 
IP=IP-1 
After every occurrence of the ATTR2 event: 
If ATTR>=QP, then schedule the ORDER() event to occur 
without delay. 
(Time ties are broken by an execution priority of 2.) 


17. The COUNT() event models the number of successful repairs 
in a batch. 
This event causes the following state change(s): 
COUNT=COUNT+ (X [0] <=RSR) 
After every occurrence of the COUNT event: 
If RQUE>0, then schedule the RPORD() event to occur without 
delay. 
If RQUE==0 and X[2]==0, then schedule the INPOS() event to 
occur without delay. 
(Time ties are broken by an execution priority of 8.) 


18. The INPOS() event models the induction of a batch of 
carcasses into repair. 
This event causes the following state change(s): 
IP=IP 
No additional events are scheduled here. 
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19. The RINIT() event models the reinitialization of TIME, L, 
T2, QP, and QR. 
This event causes the following state change(s) : 
TIME=850 
L=8.2 
T221.3 
QP=10 
QR=5 
After every occurrence of the RINIT event: 
Unconditionally, begin the next run of the simulation; 
that is, schedule the INIT() event to occur without 
delay. 
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APPENDIX F 


MODEL DEFAULTS 


Model Name: THESIS 

Model Description: DEPOT LEVEL REPAIRABLES 
Output File: THSVER8 . OUT 

Run Mode: HIGH SPEED 

Trace Vars: IP 

Random Number Seed:18645 

Initial Values: 8.2, 1.3,850, 10, 5, 10 


Ending Condition: STOP_ON_EVENT 
Event: RUN Number to Run: 1 
Trace Events: ALL EVENTS TRACED 


STATE VARIABLES 


State Variable #1 


Name: L 

Description: PROCUREMENT LEAD TIME 
Type: REAL 

Size: 1 


State Variable #2 


Name: CRR 

Description: CARCASS RETURN RATE 
Type: REAL 

Size: 1 


State Variable #3 


Name: RSR 

Description: REPAIR SURVIVAL RATE OF A CARCASS 
Type: REAL 

Size: 1 


State Variable #4 


Name: T2 

Description: REPAIR TURNAROUND TIME 
Type: REAL 

Size: 1 
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State Variable #5 


Name: D 

Description: AVERAGE QUARTERLY DEMAND 
Type: REAL 

Size: 1 


State Variable #6 


Name: QP 
Description: PROCUREMENT ORDER QUANTITY 
Type: INTEGER 
Size: 1 
State Variable #7 
Name : QR 
Description: REPAIR ORDER QUANTITY 
Type: INTEGER 
Size: 1 


State Variable #8 


Name: INV 

Description: ACTUAL ON HAND INVENTORY 
Type: INTEGER 

Size: 1 


State Variable #9 


Name: IP 

Description: INVENTORY POSITION 
Type: INTEGER 

Size: 1 


State Variable #10 


Name: QUEREP 

Description: NUMBER OF UNITS IN QUEUE FOR REPAIR 
Type: INTEGER 

Size: 1 


State Variable #11 


Name: QUEORD 

Description: NUMBER OF UNITS ON ORDER 
Type: INTEGER 

Size: 1 


State Variable #12 


Name: RQUE 

Description: NUMBER OF CARCASSES AWAITING REPAIR ORDER 
Type: INTEGER 

Size: 1 
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State Variable #13 


Name: TIME 

Description: LENGTH OF TIME FOR ONE RUN OF THE SIMULATION. 
Type: REAL 

Size: 1 


State Variable #14 


Name: Cc 
rescription: CLOCK TIME HOLDERS AND TIME WEIGHTED INVENTORY 
Type: REAL 
Size: 10 
State Variable #15 
Name : ATTR 
Description: CUMULATIVE TOTAL OF CONDEMN/LOST CARCASSES 
Type: INTEGER 
Size: 1,1 
State Variable #16 
Name: Xx 
Description: HOLDER OF A RANDOM VARIABLE FOR A STATE. 
Type: REAL 
Size: 3 


State Variable #17 


Name: REP 
Description: EVALUATION PROCESSING TIME FOR REPAIR. 
Type: REAL 
Size: 1,1 

State Variable #18 
Name: COUNT 
Description: NUMBER OF REPAIRABLE CARCASSES IN A BATCH 
Type: INTEGER 
Size: 1 

VERTICES 

Vertex #1 

Name: RUN 


Description: START OF THE SIMULATION 
State Changes: 
Parameter(s): L,T2,TIME,QP,QR 
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Vertex #2 
Name: 





DMAND 


Description: DEMAND OF A DEPOT LEVEL REPAIRABLE 
State Changes: C{6] =C [6] +C[3] *INV* (INV<0), INV=INV-1, 
X[1]=RND, IP=IP-1 
Parameter (s): 
Vertex #3 
Name: ORDER 
Description: ORDERING OF QP UNITS OF NEW MATERIAL 


State Changes: 
Parameter (s) : 


Vertex #4 
Name: 
Description: 
State Changes: 


Parameter (s) : 


Vertex #5 
Name: 
Description: 
State Changes: 


Parameter (s): 


Vertex #6 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #7 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #8 
Name: 
Description: 
State Changes: 
Parameter (s) : 


Vertex #9 
Name: 
Description: 
RUN 
State Changes: 
Parameter (s) : 


IP=IP+QP, QUEORD=QUEORD + QP 


RCVD2 

RECEIPT OF AN ORDER OF QP UNITS 

C(6] =C(6])+C[3] *INV* (INV<0), INV=INV+QP, 
QUEORD=QUEORD - QP 


CARTN 

CARCASS RETURNING TO THE STOCK POINT. 

RQUE=RQUE+1, COUNT=1, IP=IP+RQUE* (RQUE==QR), 
QUEREP=QUEREP+ (RQUE==QR) *RQUE, X[2] =0 


RPORD 
DECISION TO REPAIR OR CONDEMN A CARCASS 
X [0] =RND, QR=QR-1 


DEPOT 
DEPOT REPAIRING A CARCASS 


ATTRI 
LOSS OR CONDEMNATION OF A CARCASS 
ATTR=ATTR+1 


STATS 
GATHERING OF AVERAGES AND END OF A SIMULATION 


C[4] =C [4] /830,C[5] =C(5] /830,C[6] =C[6] /-830 
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Vertex #10 


Name: INIT 
Description: INITIALIZES IP, INV, D, CRR, RSR AND REP 
Ss t a t e Cc h a n g e 8 


INV=72, IP=INV, D=DISK{THS.DAT;0},CRR=DISK{THS.DAT;0}, 
RSR=DISK{THS.DAT;0}, REP=DISK{THS.DAT;0} 
Parameter (8s) : 


Vertex #11 
Name: ATTR1 
Description: INCREASE ATTR AND DECREASE IP BY ONE 
State Changes: QUEREP=QUEREP-1 
Parameter (s) : 


Vertex #12 
Name: RCVD1 
Description: RECORDING OF INV JUST BEFORE A RECEIPT 
State Changes: C[1]=C[2], C({2]=CLK, C[3]=(C[2])-C[1] ) * (CLK>20), 
c{[4J)=+c[4J]+c(3]*INV, 
C(5] =C(5] +C [3] * INV* (INV>=0) 
Parameter (s) : 


Vertex #13 
Name: SHIP1 
Description: RECEIPT OF A REPAIRED CARCASS 
Ss t a 3 e Cc h a n g e s 
C(6)=C[(6} +C[3] * INV* (INV<0) , QUEREP=QUEREP-1, INV=INV+1 
Parameter (s) : 


Vertex #14 
Name: DMND1 
Description: CAPTURING INV JUST PRIOR TO A DEMAND 
State Changes: C(1]=C[2), C{2]=CLK, C(3)=(C{2} -C[1} ) * (CLK>20), 
c[4J)=c[{4])+c [3] * INV, 
C(5] =C(5]+C[3] *INV* (INV>=0) 
Parameter (s) : 


Vertex #15 
Name: RPAIR 
Description: ACCUMULATION OF TIME-WEIGHTED INV VALUES 
State Changes: C[1]=C[2], C[2] =CLK, C[3]=(C[2]-C[1] ) * (CLK>20), 
cf4J)=c(4)])+cC(3]* INV, 
C[5] =C(5] +C [3] *INV* (INV>=0) 
Parameter (s) : 
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Vertex #16 
Name: ATTR2 
Description: REDUCTION OF IP DUE TO A REPAIR ATTRITION 
State Changes: ATTR=ATTR+1, IP=IP-1 
Parameter (s) : 


Vertex #17 
Name: COUNT 
Description: NUMBER OF SUCCESSFUL REPAIRS IN A BATCH 
State Changes: COUNT=COUNT+ (X [0] <=RSR) 
Parameter (s) : 


Vertex #18 
Name: INPOS 
Description: INCREASE OF IP BY INDUCTING A BATCH IN REPAIR 
State Changes: IP=IP 
Parameter (s) : 


Vertex #19 
Name: RINIT 
Description: 
State Changes: TIME=850,L=8.2,T2=1.3,QP=10,QR=5 
Parameter (s) : 


EDGES 

Edge #1 
Description: A CARCASS IS SUCCESSFULLY TURNED-IN 
Type: Scheduling 
Origin: DMAND 
Destination: CARTN 
Condition: X [1] <=CRR 
Delay: 0 
Priority: 3 
Attributes: 

Edge #2 
Description: BEGIN THE SIMULATION 
Type: Scheduling 
Origin: RUN 
Destination: INIT 
Condition: == 
Delay: 0 
Priority: 5 


Attributes: 








Edge #3 


Description: SCHEDULE THE END OF THE RUN 
Type: Scheduling 
Origin: INIT 
Destination: STATS 
Condition: l==1 
Delay: TIME 
Priority: 4 
Attributes: 

Edge #4 
Description: A CARCASS IS NOT RETURNED TO THE STOCK POINT 
Type: Scheduling 
Origin: DMAND 
Destination: ATTRI 
Condition: CRR<X [1] 
Delay: 0 
Priority: 3 
Attributes: 

Edge #5 
Description: PROCURE QP UNITS 
Type: Scheduling 
Origin: ATTRI 
Destination: ORDER 
Condition: ATTR>=QP 
Delay: 0 
Priority: 3 
Attributes: 

Edge #6 
Description: SCHEDULE THE RECEIPT OF QP UNITS 
Type: Scheduling 
Origin: ORDER 
Destination: RCVD1 
Condition: l== 
Delay: L 
Priority: 5 
Attributes: 

Edge #7 
Description: ADD QP UNITS TO INVENTORY (INV) 
Type: Scheduling 
Origin: RCVD1 
Destination: RCVD2 
Condition: == 
Delay: * 
Priority: 5 
Attributes: 


124 


Edge #8 


Description: SCHEDULE ANOTHER DEMAND 
Type: Scheduling 
Origin: DMAND 
Destination: DMND1 
Condition: 1==1 
Delay: - (1/D) *LN{RND} 
Priority: 5 
Attributes: 
Edge #9 
Description: RECORD INV AND IP JUST PRIOR TO A DEMAND 
Type: Scheduling 
Origin: DMND1 
Destination: DMAND 
Condition: a= 
Delay: * 
Priority: 5 
Attributes: 
Edge #10 
Description: SHIP THE RFI CARCASS TO THE STOCK POINT 
Type: Scheduling 
Origin: DEPOT 
Destination: RPAIR 
Condition: l== 
Delay: T2 
Priority: 3 
Attributes: 
Edge #11 
Description: SCHEDULE THE REPAIR ATTRITION 
Type: Scheduling 
Origin: RPORD 
Destination: ATTR1 
Condition: X[0] >RSR 
Delay: REP* COUNT 
Priority: 3 
Attributes: 
Edge #12 
Description: SCHEDULE THE REPAIR OF A NRFI UNIT 
Type: Scheduling 
Origin: RPORD 
Destination: DEPOT 
Condition: X[0]<=RSR 
Delay: REP* COUNT 
Priority: a 
Attributes: 
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Edge #13 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 

Ar ributes: 


Edge #14 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #15 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #16 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 


Edge #17 
Description: 
Type: 

Origin: 
Destination: 
Condition: 
Delay: 
Priority: 
Attributes: 





SEND A BATCH OF QP NRFI CARCASSES TO DEPOT 
Scheduling 
CARTN 
RPORD 
RQUE>=QR 
0 
5 


INCREASE ATTR AND REDUCE IP BY ONE 
Scheduling 
ATTR1 
ATTR2 
1== 
0 
7 


INCREASE THE COUNT OF REPAIRABLES IN A BATCH 
Scheduling 
RPORD 
COUNT 
1l==1 
0 
5 


SCHEDULE THE REVIEW OF THE NEXT CARCASS 
Scheduling 

COUNT 
RPORD 


PROCURE QP UNITS 
Scheduling 
ATTR2 
ORDER 
ATTR>=QP 
0 
2 
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Edge #18 


Description: RECEIVE A REPAIRED UNIT 
Type: Scheduling 
Origin: RPAIR 
Destination: SHIP1 
Condition: l==1 
Delay: * 
Priority: 5 
Attributes: 

Edge #19 
Description: BEGIN THE NEXT RUN OF THE SIMULATION. 
Type: Scheduling 
Origin: STATS 
Destination: RINIT 
Condition: ==1&SET{18645} 
Delay: 0 
Priority: 8 
Attributes: 

Edge #20 
Description: BEGIN THE NEXT RUN OF THE SIMULATION. 
Type: Scheduling 
Origin: RINIT 
Destination: INIT 
Condition: l==1 
Delay: 0 
Priority: 5 
Attributes: 

Edge #21 
Description: BEGIN THE DEMAND CYCLE OF EACH RUN 
Type: Scheduling 
Origin: INIT 
Destination: DMND1 
Condition: 1==1 
Delay: 0 
Priority: 5 
Attributes: 

Edge #22 
Description: 
Type: Scheduling 
Origin: CARTN 
Destination: INPOS 
Condition: RQUE==QOR 
Delay: 0 
Priority: 8 
Attributes: 
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APPENDIX G 


The SIGMA Model, THESIS, is a discrete event simulation. It 
models depot level repairables. 


I. STATE VARIABLE DEFINITIONS. 


For this simulation, the following state variables are 
defined: 
L: procurement lead time (real valued) 
CRR: carcass return rate (real valued) 
RSR: repair survival rate of a carcass. (real valued) 
T2: repair turnaround time (real valued) 
D: average quarterly demand (real valued) 


QP: procurement order quantity (integer valued) 

QR: repair order quantity (integer valued) 

INV: actual on hand inventory (integer valued) 

IP: inventory position (integer valued) 

QUEREP: number of units in queue for repair (integer 
valued) 

QUEORD: number of units on order (integer valued) 


RQUE: number of carcasses awaiting repair order 
(integer valued) 
TIME: length of time for one run of the simulation. 
(real valued) 
C{10]: clock time holders and time weighted inventory 
(real valued) 
ATTR: cumulative total of condemn/lost carcasses 
(integer valued) 
X(3]: holder of a random variable for a state. (real 
valued) 
REP: evaluation processing time for repair. (real valued) 
COUNT: number of repairable carcasses in a batch 
(integer valued) 


II. EVENT DEFINITIONS. 


Simulation state changes are represented by event vertices 
(nodes or balls) in a SIGMA graph. Event vertex parameters, 
if any, are given in parentheses. 

Logical and dynamic relationships between pairs of events 
are represented in a SIGMA graph by edges (arrows) between 
event vertices. Unless otherwise stated, vertex execution 
priorities, to break time ties, are equal to 5. 
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1. The RUN(L,T2,TIME,QP,QR) event models the start of the 
simulation. 

Initial values for, L,T2,TIME,QP,QR, are needed for each 
run. 
After every occurrence of the RUN event: 
Unconditionally, begin the simulation; 

that is, schedule the INIT() event to occur without 

delay. 


2. The DMAND() event models the demand of a depot level 
repairable. 

This event causes the following state change(s) : 

C(6] =C(6] +C[3) * INV* (INV<0) 
INV=INV-1 

X [1] =RND 

IP=IP-1 

After every occurrence of the DMAND event: 

If X[1])<=CRR, then a carcass is successfully turned-in; 
that is, schedule the CARTN() event to occur without 
delay. 

(Time ties are broken by an execution priority of 3.) 

If CRR<X[1], then a carcass is not returned to the stock 

point; 
that is, schedule the ATTRI() event to occur without 
delay. 
(Time ties are broken by an execution priority of 3.) 

Unconditionally, schedule another demand; 
that is, schedule the DMND1() event to occur in 
-(1/D)*LN{RND} time units. 


3. The ORDER() event models the ordering of QP units of new 
material. 
This event causes the following state change(s) : 
IP=IP+QP 
QUEORD=QUEORD + QR 
After every occurrence of the ORDER event: 
Unconditionally, schedule the receipt of QP units; 
that is, schedule the RCVD1() event to occur in L 
time units. 





4. The RCVD2() event models the receipt of an order of QP 
units. 
This event causes the following state change(s) : 
C(6) =C(6] +C[3] *INV* (INV<0) 
INV=INV+QP 
QUEORD=QUEORD-QP 
No additional events are scheduled here. 
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5. The CARTN() event models the carcass returning to the stock 
point. 
This event causes the following state change(s): 
RQUE=RQUE+1 
COUNT=1 
IP#IP+RQUE* (RQUE==QR) 
QUEREP=QUEREP+ (RQUE==QR) *RQUE 
X(2] =0 
After every occurrence of the CARTN event: 
If RQUE>=QR, then send a batch of QP NRFI carcasses to 
depot; 
that is, schedule the RPORD() event to occur without 
delay. 
If RQUE==QOR, then schedule the INPOS() event to occur 
without delay. 
(Time ties are broken by an execution priority of 8.) 


6. The RPORD() event models the decision to repair or condemn 
a carcass. 
This event causes the following state change(s) : 
X [0] =RND 
QR=QR-1 
After every occurrence of the RPORD event: 
If X[0]>RSR, then schedule the repair attrition; 
that is, schedule the ATTR1() event to occur in 
REP*COUNT time units. 
(Time ties are broken by an execution priority of 3.) 
If X[0]<=RSR, then schedule the repair of a NRFI unit; 
that is, schedule the DEPOT() event to occur in 
REP*COUNT time units. 
(Time ties are broken by an execution priority of 3.) 
Unconditionally, increase the count of repairables in a 
batch; 
that is, schedule the COUNT() event to occur without 
delay. 


7. The DEPOT() event models the depot repairing a carcass. 
After every occurrence of the DEPOT event: 
Unconditionally, ship the RFI carcass to the stock point; 

that is, schedule the RPAIR() event to occur in T2 
time units. 
(Time ties are broken by an execution priority of 3.) 
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8. The ATTRI() event models the loss or condemnation of a 
carcass. 
This event causes the following state change(s): 
ATTRsATTR+1 
After every occurrence of the ATTRI event: 
If ATTR>=QP, then procure QP units; 
that is, schedule the ORDER() event to occur without 


delay. 
(Time ties are broken by an execution priority of 3.) 


9. The STATS() event models the gathering of averages and end 
of a simulation run. 
This event causes the following state change(s): 
C[4) =C [4] /830 
C(5) =C (5) /830 
C{6} =C(6) /-830 
After every occurrence of the STATS event: 
If 1l=*=1 and SET{18645}, then begin the next run of the 
simulation; that is, schedule the RINIT() event to occur 
without delay. 
(Time ties are broken by an execution priority of 8.) 


10. The INIT() event models the initializes IP, INV, D, CRR, 
RSR, and REP. 
This event causes the following state change(s): 
INV=72 
IP=INV 
D=DISK{THS.DAT;0} 
CRR=DISK{ THS.DAT;0 
RSR=DISK{ THS . DAT; 0 
REP=DISK{THS .DAT; 0 
After every occurrence of the INIT event: 
Unconditionally, schedule the end of the run; 
that is, schedule the STATS() event to occur in TIME 
time units. 
(Time ties are broken by an execution priority of 4.) 
Unconditionally, begin the demand cycle of each run; 
that is, schedule the DMND1() event to occur without 
delay. 








11. The ATTR1() event models the increase ATTR and decrease IP 
by one. 
This event causes the following state change(s): 
QUEREP=QUEREP-1 
After every occurrence of the ATTR1 event: 
Unconditionally, increase attr and reduce IP by one; 
that is, schedule the ATTR2() event to occur without 
delay. 
(Time ties are broken by an execution priority of 7.) 
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12. The RCVD1() event models the recording of INV just before 
a receipt. 
This event causes the following state change(s): 
C(1}) =C[2) 
C(2) =CLK 
C(3} =(C(2) -C[1) ) * (CLK>20) 
C{4) =C [4] +C[3] * INV 
C(S] =C(5) +C[3] *INV* (INV>=0) 
After every occurrence of the RCVD1 event: 
Unconditionally, add QP units to inventory (INV); 
that is, immediately execute the RCVD2() event. 


13. The SHIP1() event models the receipt of a repaired 
carcass. 
This event causes the following state change(s): 
C[{6) =C (6) +C[3] *INV* (INV<0) 
QUEREP=QUEREP-1 
INV=INV+1 
No additional events are scheduled here. 


14. The DMND1() event models the capturing INV just prior to 
a demand. ‘ 
This event causes the following state change(s) : 
C{1) =C[2] 
C(2] =CLK 
C{3} =(C[2] -C[1] ) * (CLK>20) 
C [4] =C(4]+C(3] * INV 
C[5] =C([5]) +C (3) *INV* (INV>=0) 
After every occurrence of the DMND1 event: 
Unconditionally, record inv and ip just prior to a demand; 
that is, immediately execute the DMAND() event. 


15. The RPAIR() event models the accumulation of time-weighted 
INV values. 
This event causes the following state change(s) : 
C[1] =c[2) 
C[2) =CLK 
C[3) =(C[2] -C[1) ) * (CLK>20) 
C(4] =«C [4] +C[3] * INV 
C[(5} =C(5] +C [3] *INV* (INV>=0) 
After every occurrence of the RPAIR event: 
Unconditionally, receive a repaired unit; 
that is, immediately execute the SHIP1() event. 


16. The ATTR2() event models the reduction of IP due to a 
repair attrition. 
This event causes the following state change(s) : 
ATTRzATTR+1 
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IPsIP-1 
After every occurrence of the ATTR2 event: 
If ATTR>=QP, then procure QP units; 
that is, schedule the ORDER() event to occur without 
delay. 
(Time ties are broken by an execution priority of 2.) 


17. The COUNT() event models the number of successful repairs 
in a batch. 
This event causes the following state change(s) : 
COUNT=COUNT+ (X [0] <=RSR) 
After every occurrence of the COUNT event: 
If RQUE>0, then schedule the review of the next carcass; 
that is, schedule the RPORD() event to occur without 
delay. 


18. The INPOS() event models the increase of IP by inducting 
a batch in repair. 
This event causes the following state change(s) : 
IP=#IP 
No additional events are scheduled here. 


19. The RINIT() event: 

This event causes the following state change(s): 
TIME=850 
L=8.2 
T22=1.3 
QP=10 
QR=5 

After every occurrence of the RINIT event: 

Unconditionally, begin the next run of the simulation.; 
that is, schedule the INIT() event to occur without 
delay. 
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woe we wmew ewww eer eo cee ewer nwo ewer er een errr renew 


CRR RSR 
0.00 1.00 
0.20 0.00 
0.20 0.00 
0.20 0.00 
0.20 0.00 
0.20 0.00 
0.20 0.00 
0.20 0.20 
0.20 0.20 
0.20 0.20 
0.20 0.20 
0.20 0.20 
0.20 0.20 
0.20 0.40 
0.20 0.40 
0.20 0.40 
0.20 0.40 
0.20 0.40 
0.29 0.40 
0.20 0.60 
0.20 0.60 
0.20 0.60 
0.20 0.60 
0.20 0.60 
0.20 0.60 
0.20 0.80 
0.20 0.80 
0.20 0.80 
0.20 0.80 
0.20 0.80 
0.20 0.80 
0.20 1.00 
0.20 1.00 
0.20 1.00 
0.20 1.00 
0.20 1.00 
0.20 1.00 
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518.364 


1.384 
1.384 
1.384 
1.384 
1.384 
1.384 


1.382 
1.699 
1.719 
2.458 
2.820 
2.157 


2.681 
3.195 
3.529 
3.136 
3.874 
4.068 


1.013 
1.275 
2.117 
2.330 
2.321 
2.549 


2.304 
1.987 
2.551 
1.911 
2.370 
1.756 


2.972 
2.972 
2.972 
2.972 
2.972 
2.972 


1.0000 


0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


0.0000 
0.0001 
0.0001 
0.0007 
0.0015 
0.0003 


0.0011 
0.0029 
0.0048 
0.0026 
0.0075 
0.0096 


0.0000 
0.0000 
0.0003 
0.0005 
0.0005 
0.0009 


0.0005 
0.0002 
0.0009 
0.0002 
0.0006 
0.0001 


0.0020 
0.0020 
0.0020 
0.0020 
0.0020 
0.0020 





wwreeewewr nw erewrecenwn2nnewerzrrznrzrececewnenrenenwere2ezrerecerer eee cece 


0.40 0.00 0.00 1.699 0.0001 
0.40 0.00 0.25 1.699 0.0001 
0.40 0.00 0.50 1.699 0.0001 
0.40 0.00 0.75 1.699 0.0001 
0.40 0.00 1.00 1.699 0.0001 
0.40 0.00 1.30 1.699 0.0001 
0.40 0.20 0.00 2.254 0.0004 
0.40 0.20 0.25 1.620 0.0001 
0.40 0.20 0.50 1.697 0.0001 
0.40 0.20 0.75 2.097 0.0003 
0.40 0.20 1.00 1.744 0.0001 
0.40 0.20 1.30 1.478 0.0000 
0.40 0.40 0.00 2.672 0.0011 
0.40 0.40 0.25 3.693 0.0060 
0.40 0.40 0.50 3.557 0.0049 
0.40 0.40 0.75 3.841 0.0072 
0.40 0.40 1.00 4.160 0.0106 
0.40 0.40 1.30 3.325 0.0035 
0.40 0.60 0.00 1.270 0.0000 
0.40 0.60 0.25 0.643 0.0000 
0.40 0.60 0.50 1.024 0.0000 
0.40 0.60 0.75 1.206 0.0000 
0.40 0.60 1.00 1.151 0.0000 
0.40 0.60 1.30 1.216 0.0000 
0.40 0.80 0.00 1.279 0.0000 
0.40 0.80 0.25 1.090 0.0000 
0.40 0.80 0.50 1.200 0.0000 
0.40 0.80 0.75 0.579 0.0000 
0.40 0.80 1.00 0.919 0.0000 
0.40 0.80 1.30 1.349 0.0000 
0.40 1.00 0.00 1.517 0.0000 
0.40 1.00 0.25 1.517 0.0000 
0.40 1.00 0.50 1.517 0.0000 
0.40 1.00 0.75 1.517 0.0000 
0.40 1.00 1.00 1.517 0.0000 
0.40 1.00 1.30 1.517 0.0000 
0.60 0.00 0.00 1.098 0.0000 
0.60 0.00 0.25 1.098 0.0000 
0.60 0.00 0.50 1.098 0.0000 
0.60 0.00 0.75 1.098 0.0000 
0.60 0.00 1.00 1.098 0.0000 
0.60 0.00 1.30 1.098 0.0000 


142 








143 


1. 
2 
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0.80 0.40 0.00 1.431 0.0000 
0.80 0.40 0.25 1.032 0.0000 
0.80 0.40 0.50 1.202 0.0000 
0.80 0.40 0.75 1.320 0.0000 
0.80 0.40 1.00 1.427 0.0000 
0.80 0.40 1.30 2.484 0.0008 
0.80 0.60 0.00 3.350 0.0037 
0.80 0.60 0.25 1.741 0.0001 
0.80 0.60 0.50 0.929 0.0000 
0.80 0.60 0.75 1.297 0.0000 
0.80 0.60 1.00 1.541 0.0001 
0.80 0.60 1.30 1.869 0.0002 
0.80 0.80 0.00 0.775 0.0000 
0.80 0.80 0.25 0.957 0.0000 
0.80 0.80 0.50 0.855 0.0000 
0.80 0.80 0.75 0.913 0.0000 
0.80 0.80 1.00 1.218 0.0000 
0.80 0.80 130 2.380 0.0006 
0.80 1.00 0.00 2.321 0.0005 
0.80 1.00 0.25 2.321 0.0005 
0.80 1.00 0.50 2.321 0.0005 
0.80 1.00 0.75 2.321 0.0005 
0.80 1.00 1.00 2.321 0.0005 
0.80 1.00 1.30 2.321 0.0005 
1.00 0.00 0.00 540.996 1.0000 
1.00 0.00 0.25 540.996 1.0000 
1.00 0.00 0.50 540.996 1.0000 
1.00 0.00 0.75 540.996 1.0000 
1.00 0.00 1.00 540.996 1.0000 
1.00 0.00 1.30 540.996 1.0000 
1.00 0.20 0.00 2.358 0.0006 
1.00 0.20 0.25 0.662 0.0000 
1.00 0.20 0.50 1.350 0.0000 
1.00 0.20 0.75 0.983 0.0000 
1.00 0.20 1.00 0.498 0.0000 
1.00 0.20 1.30 1.555 0.0001 
1.00 0.40 0.00 0.700 0.0000 
1.00 0.40 0.25 0.984 0.0000 
1.00 0.40 0.50 1.460 0.0000 
1.00 0.40 0.75 1.421 0.0000 
1.00 0.40 1.00 1.455 0.0000 
1.00 0.40 1.30 1.656 — 0.0001 
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APPENDIX J 


Time Event Count CRR RSR REP SAFTY 
850.000 STATS 1 0 i 1.3 -12.4 
850.000 STATS 1 0.2 0 0 -15.1 
850.000 STATS 1 0.2 0 0.25 -15.1 
850.000 STATS 1 0.2 0 0.5 -15.1 
850.000 STATS 1 0.2 0 0.75 -15.1 
850.000 STATS 1 0.2 0 1 -15.1 
850.000 STATS 1 0.2 0 1.3 “15.1 
850.000 STATS 1 0.2 0.2 0 -11.7 
850.000 STATS 1 0.2 0.2 0.25 -11.8 
850.000 STATS 1 0.2 0.2 0.5 e119 
850.000 STATS 1 0.2 0.2 0.75 “12.1 
850.000 STATS 1 0.2 0.2 1 -12.2 
850.000 STATS 1 0.2 0.2 ba -12.3 
850.000 STATS 2 0.2 0.4 0 -8.82 
850.000 STATS 1 0.2 0.4 0.25 -8.74 
850.000 STATS 1 0.2 0.4 0.5 -8.98 
850.000 STATS 1 0.2 0.4 0.75 -9.25 
850.000 STATS 1 0.2 0.4 1 -9.59 
850.000 STATS 1 0.2 0.4 i.3 -9.98 
850.000 STATS i 0.2 0.6 0 -6.36 
850.000 STATS 1 0.2 0.6 0.25 -6.09 
850.000 STATS 1 0.2 0.6 0.5 -6.52 
850.000 STATS a 0.2 0.6 0.75 -7.08 
850.000 STATS 1 0.2 0.6 1 -7.59 
850.000 STATS zt 0.2 0.6 1.3 -8.2 
850.000 STATS 1 0.2 0.8 0 -3.78 
850.000 STATS 1 0.2 0.8 0.25 -3.34 
850.000 STATS 1 0.2 0.8 0.5 -4.02 
850.000 STATS 1 0.2 0.8 0.75 -4.73 
850.000 STATS a 0.2 0.8 1 -5.53 
850.000 STATS 1 0.2 0.8 1.3 -6.2 
850.000 -.ATS 1 0.2 1 0 -1 Gh 
850.000 STATS 1 0.2 1 0.25 -0.523 
850.000 STATS 1 0.2 1 0.5 -1.33 
850.000 STATS 1 0.2 1 0.75 -2.3 
850.000 STATS 1 0.2 1 1 -3.14 
850.000 STATS 1 0.2 1 1.3 -4.19 
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650.000 
850.000 
850.000 
850.000 
850.000 
850.000 


850.000 
850.000 
850.000 
850.000 
850.000 
850.000 


850.000 
850.000 
850.000 
850.000 
850.000 
850.000 


850.000 
850.000 
850.000 
850.000 
850.000 
850.000 


850.000 
850.000 
850.000 
850.000 
850.000 
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